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1. Introduction 

The last decade can be characterized by an impressive diversity of techniques mon- 
itoring the artificial and natural satellite dynamics, as well as the Earth rotation: im- 
proved laser technology, renewed Doppler techniques, satellite radar altimetry, massive 
usage of the Global Positioning System (GPS), etc. Each of these techniques is opti- 
mally tailored to a specific type of application or scientific problem. For example, it 
appears that laser tracking (SLR: see on the WEB: http:/ /ilrs. gsfc.nasa.gov) of pas- 
sive geodynamics satellites (LAGEOS, LAGEOS II, Starlette, Stella, Ajisai, Etalon I 
and II) over relatively long time intervals provides an excellent method for determin- 
ing the long-term variations of the geopotential [1] (including tidal effects) and many 
small non-gravitational phenomena [2]. It turns out that the precision reached in the 
latest years by such technique in measuring the position of the passive laser-ranged 
geodetic satellites LAGEOS and LAGEOS II amounts to 1 cm or better(^). See, e.g., 
http: / / earth.agu.org/revgeophys / marshaOl/ nodel.html 

Such astonishing levels of accuracy have disclosed an unexpectedly wide field of re- 
search in space geodesy, geophysics and fundamental physics. E.g., now it is possible 
to plan satellite-based experiments devoted to the experimental control even of some 
tiny post-Newtonian features of the Earth'gravitational field predicted by Einstein' s 
General Relativity [3]. They are usually expressed in terms of certain solve-for least 
squares fits' parameters. In evaluating the precision of the results of these experiments 
it must be considered that, in general, the main error does not consist of the standard 
statistical error of the fits but of the various systematic errors. They are induced by a 
complete set of other physical effects acting upon the satellites to be employed [4] . Such 
perturbations are often quite larger than the relativistic effects investigated and induce 
systematic errors to be correctly evaluated and assessed. Many of these aliasing effects 
are traditionally investigated by geophysics and space geodesy, so that such experiments 
are multidisciplinary efforts which cover many scientific areas until now separated [5, 6]. 
Indeed, the expansion of the SLR network, together with improved system accuracies, 
has enabled the laser data to contribute directly to improving orbit force models. It 
allows one to calculate more accurately, among other things, the error budget of many 
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space-based general relativistic experiments: this is the main topic of the present work. It 
does not treat the measurement modelling errors specifically related to the laser-ranging 
technique. 

I'l. The Lense-Thirring drag of inertial frames. - One of the most interesting topic 
in General Relativity is the structure of the spacetime around a spherically symmetric 
rotating mass-energy distribution. Indeed, in the slow-motion and weak-field approxima- 
tion, it exhibits the characteristic feature of exerting a non-central force on a test particle 
due to the total angular monientiini of the central object, contrary to Newtonian me- 
chanics in which the gravitational action of a body is caused only by its mass, regardless 
to its rotational motion. Because of the formal analogies with the electrodynamics, this 
effect, deduced from the equations of Einstein for the first time by Lense and Thirring in 
1918 [7], is also defined as gravitomagnetism. A comprehensive review of its properties 
can be found in [8, 9, 10]. 

In the past few years we have seen increasing efforts, both from a theoretical and an 
experimental point of view, devoted to the measurement of the Lense-Thirring effect in 
the weak gravitational field of the Earth by means of artificial satellites. At present, 
there arc two main proposals which point towards the implementation of this goal: the 
Gravity Probe-B mission [11], and the approach proposed in [12] which consists in using 
the actually orbiting LAGEOS laser-ranged satellite and launching another satellite of 
LAGEOS type, the LARES, with the same orbital parameters of LAGEOS except for 
the inclination, which should be supplementary with respect to it, the eccentricity which 
should be one order of magnitude larger, and the area-to-mass ratio which should be 
smaller so to reduce the impact of the non-gravitational orbital perturbations (see Ap- 
pendix A for the LAGEOS type satellites' data). The observable would be a secular linear 
trend built up with the sum of the residuals of the longitudes of the ascending nodes fl 
of LARES and LAGEOS (see section 2 and fig. 1 for an explanation of the Keplerian 
orbital elements in terms of which the orbit of a satellite can be parameterized). The 
proposed orbital geometry would allow one to minimize the impact of the aliasing trend 
due to the mismodelling in the classical nodal precessions generated by the oblateness of 
the Earth, which would represent the main systematic error. 

The GP-B mission is aimed to the detection of the motion of a spinning particle. At 
present, both the GP-B and the LARES have not yet been launched: however, while the 
GP-B is scheduled to fiy in fall 2002, the fate of LARES project is still uncertain. 

Recently, Ciufolini, in ref.[14], has put forward an alternative strategy based on the 
utilization of the already existing LAGEOS and LAGEOS II which allowed the detection 
of the Lense-Thirring drag at a precision level of the order of 20% [15, 16]. While the 
GP-B mission is focused on the gravitomagnetic dragging of the spin of a freely falling 
body, in the LAGEOS experiment it is the entire orbit of the satellite which is considered 
to undergo the secular Lense-Thirring precession. More exactly, among the Keplerian 
orbital elements [17, 18] , the node fl and the perigee u are affected by the gravitomagnetic 
perturbation. For LAGEOS and LAGEOS II it amounts to 



(1) fjLAGEOS ^ 3^ j^^g/y^ 

(2) f^LAGEOS II ^ 31 5 ^^g/y^ 

(3) wj^^^'^os ~ 31.6 mas/y, 

(4) a;L^GEOS ii ^ ^^^/^^ 
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where mas/y stands for milliarcseconds per year. 

In section 2 a brief derivation of such results is presented. It fohows ordinary Unear 

satellite perturbation theory. Indeed, in the slow-motion and weak-field approximation, 
the gravitomagnetic potential is treated as a classical disturbing term with respect to 
the Newtonian gravitoelectric monopole. 

The observable quantity proposed in [14] is a suitable combination of the orbital 
residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II 

(5) y = Jfi^p + ci X Stl^^p + C2 X Sojgp ~ /xlt x 60.2. 

In it Ci = 0.295, C2 = —0.35, /Ult is the scaling, solve-for parameter, equal to 1 in General 
Relativity and in classical mechanics, to be determined and (5rJj,j,p, <5f^oxp' '^^oxp 
the orbital residuals, in mas, calculated with the aid of some orbit determination software 
like UTOPIA or GEODYN, of the nodes of LAGEOS and LAGEOS II and the perigee 
of LAGEOS II. The residuals account for any unmodelled or mismodelled physical phe- 
nomena acting on the observable analyzed. By dealing with the gravitomagnetism as an 
unmodelled physical effect. General Relativity predicts a linear trend with a slope of 60.2 
mas/y 

(6) yLT = (31 mas/y) + ci x (31.5 mas/y) + C2 x (-57 mas/y) ~ 60.2 mas/y. 

The coefficients ci and C2 have been obtained by solving for the three unknowns 5J2, 
SJi and /Ult a nonhomogeneous algebraic linear system of three equations expressing the 
observed mismodelled classical precessions of the nodes of LAGEOS and LAGEOS II 
and the perigee of LAGEOS II (see ref.[14]). They depend on the orbital parameters of 
LAGEOS and LAGEOS II and are built up so to cancel out all the static and dynamical 
contributions of degree I = 2,4 and order m = of the Earth' s gravitational field. This 
cancellation is required to reduce especially the impact of the first two even (Z = 2,4) 
zonal (m = 0) harmonics of the static geopotential. Indeed, at the present level of 
accuracy of the terrestrial gravitational field [19], the mismodelled parts of the classical 
orbital precessions of the node and the perigee, caused by the Earth's oblateness and 
parameterized with the first two even zonal coefficients J2 and J4 [1], for a single satellite 
amount to a significative part of the corresponding Lense-Thirring effect itself on the 
considered orbital element. The combined residuals of eq. (5) are affected by a complete 
set of gravitational and non-gravitational perturbations; among the former ones those 
generated by the solid Earth and ocean tides [20, 21] play an important role. 

In section 3 the amplitudes and the periods of the most relevant tidal perturbations 
acting upon the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II have 
been worked out [22]. This analysis covers both the solid I = 2 and the ocean Z = 2, 3, 4 
part of the tidal spectrum, includes the first order contributions, in the sense of the 
Kaula' s theory of perturbations [1], and also a calculation of the mismodelling in the 
first order perturbative amplitudes compared to the Lense-Thirring drag over 4 years. 

Section 4 is devoted to the study of the impact of the solid Earth and ocean tidal 
perturbations on eq. (5) with particular care to the systematic errors induced by them 
on the measurement of /zlt [23, 24]. The left hand side of eq. (5) has been calculated 
by adopting the nominal values of the perturbative amplitudes worked out in section 
3 in order to test if the first two even degree zonal tidal lines do not really affect the 
combined residuals. This is particularly important for the 18.6-year tidal constituent 
whose period is that of the longitude of the ascending node of the Moon, i.e. 18.6 years. 



4 



L. lORIO 



Indeed, over observational time spans of few years it could resemble a trend as well 
and its mismodelled part, accounted for by the residuals combination, could corrupt the 
measurement of /ilt- The rcsiilts obtained in section 3 have been employed also in order 
to obtain a simulated residual curve on which several tests have been performed for the 
various tesseral (m = 1) and sectorial (m = 2) tidal constituents, not cancelled out by 
the combined residuals, with the aim of assessing their contribution on the systematic 
uncertainty in the measurement of /ilt- 

1'2. The gravitoelectric pericenter shift. - In section 5 the experience gained in the 
Lense-Thirring LAGEOS experiment is used to propose a possible experiment for mea- 
suring the general relativistic pericenter shift of a test body due to the Schwarzschild's 
gravitoelectric part of the terrestrial gravitational field by means of the SLR data to 
LAGEOS and LAGEOS II. The approach is similar to that of the current Lense-Thirring 
LAGEOS experiment and it is based on another suitable combination of orbital residuals 
of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The possibil- 
ity of adopting different observables is examined as well. A discussion on the obtainable 
values for the Parameterized Post Newtonian parameters /3 and 7 is also included. 

1'3. The LARES mission. In section 6 we reanalyze the original LARES mission 
from the point of view of the sensitivity of the related observable to the departures of the 
LARES orbital parameters from their nominal values due to the orbital injection errors. 
Moreover, we propose an alternative combination of orbital residuals involving also the 
LAGEOS II which would yield to a smaller and more stable value of the error due to the 
mismodelling in the even zonal harmonics of the geopotential which turns out to be the 
main source of systematic error. 

1'4. Conclusions. - Section 7 is devoted to the conclusions and to the recommenda- 
tions for future work. 

2. — The Lense-Thirring effect on the orbit of a test body 

2'1. Introduction. - In this section an alternative strategy is presented in order to 

derive the Lense-Thirring effect [25]. It reveals itself useful especially in the prediction of 
the behavior of all the Keplerian orbital elements of the test body in the gravitational field 
of a rotating source. Our calculations are based on the Lagrangian planetary equations 
[1, 26, 27] and a non-central hamiltonian term whose existence, if from one hand can be 
rigorously deduced, from the other hand can be intuitively guessed by analogy from the 
corresponding term in electrodynamics for the Lagrangian of a charged particle acted 
upon by electric and magnetic fields. 

For a perfectly spherical source the well known Lense-Thirring rate equations for the 
secular precessions of the node and the perigee are reobtained. 

2'2. The gravitomagnetic potential. - In general, it can be proved [9] that the general 

relativistic equations of motion of a test particle of mass m freely falling in a stationary 
gravitational field, in the weak-field and slow-motion approximation, are given by 



In eq. (7) Eg and Bg are the gravitoelectric and gravitomagnetic fields, respectively. If a 
perfectly spherically symmetric rotating body is assumed as gravitational source, in eq. 



(7) 
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(7) Eg = —GMir/r"^ + 0{c ^) is the Newtonian gravitational field of a spherical body, 
with M its mass and G the Newtonian gravitational constant(^), while Bg is given by 
[9] 



(8) 



Bg = V X Ag 



G 



J — 3(J • ir)ir 



in which 



• / N „G J X r 
(9) Ag(r) ^ -2-—^. 

In eq. (9) J is the angular momentum of the central body. The field Ag = (/iqi , /io2j ^03)) 
named gravitomagnetic potential, is due to the off-diagonal components of the spacetime 
metric tensor 



(10) gij.u = Vniy + = 0, 1, 2, 3, 

where rj^i, is the Minkowski metric tensor. In obtaining eq. (9) a non-rotating reference 
frame K{x, y, z} with the z axis directed along the intrinsic angular momentum of the 
central spinning body J and the {x, y} plane coinciding with its equatorial plane has 
been used. The origin is located at the center of mass of the central body. Eq. (7) holds 
if Ag is not time-varying. 

The gravitomagnetic potential generates a non-central gravitational contribute due 
uniquely to the angular momentum of the gravitational source that the Newtonian me- 
chanics does not predict, though the conditions of validity of eq. (7) are the same for 
which the latter holds as well.(^) So it is possible to speak of mass-energy currents whose 
motion exerts a non-central gravitational force on a test massive body analogous to the 
Lorentz force felt by a charged particle when it is moving in a magnetic field. Indeed, its 
equations of motions 

(11) m-^.,(E+-xB) 

are formally analogous to eq. (7). Eq. (11) can be derived by means of the Lagrangian 



(12) Ce.m. = lmv''-qV + ^vA). 

2 c 

where v is the velocity of the particle and V is its scalar potential. 



(^) The terms of order 0{c^'^) are the post-Newtonian corrections due to the static 
Schwarzschild part of the metric. 

(^) Incidentally, it may be interesting to notice that eq. (7) are consistent with the fundamental 
Einstein assumption [28] that a non-accelerated reference frame with a gravitational field is 
equivalent, within certain limits, to an accelerated one without any gravitational field. Indeed, 
if a reference frame solidal with the rotating body is assumed, the equations of motion for a 
test particle arc formally identical to eq. (7): the gravitomagnetic force term ^(v x Bg) is 
substituted by the Coriolis force term 2m (v x fJ) where J7 is the angular velocity vector of the 
rotating body [29]. See also [30]. 
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Since one of the most promising ways to detect the gravitomagnetic precession consists 
of employing artificial Earth satellites, it would be helpful to derive the rate equations 
for the change in the parameters that characterize their orbits. To this aim one could 
introduce "by hand" a perturbative term k (v • Ag) in the gravitational Lagrangian of 
the particle and use it in some particular perturbative scheme; the constant k would 
be determined by means of dimensional considerations and taking into account that it 
should be built up of universal constants. In fact it is possible to show that a non-central 
term analogous to ' (v • A) can be rigorously deduced in the Lagrangian of a test body in 
the gravitational field of a spinning mass-energy distribution, and that it can be exploited 
in deriving straightforwardly the effect of the gravitomagnetic potential on the Keplerian 
orbital elements of the test body. 

2"3. The rate equations for the Keplerian orbital elements. - The relativistic La- 
grangian for a free particle in a gravitational field can be cast into the form: 

(13) £ = £(°)+£W. 

In eq. (13) the term £^^^ is to be intended as containing the contributions of the off- 
diagonal terms of the metric 

(14) £(1) = "^gokx'^xK 

In this case, recalling that the slow motion approximation is used, the eq. (14) becomes 

(15) £W=£g„ = ^(Ag.v). 

For different derivations of Cgm see also [31, 32, 26]. In this section it is proposed to 
adopt Cgm in order to deriving the Lcnsc-Thirring effect on the orbital elements of a 
particle in the field of a rotating gravitational source. 

To this aim it must be assumed that under the gravitomagnetic force the departures 
of the test body' s trajectory from the unperturbed Keplerian ellipse arc very small in 
time. This allows to introduce the concept of osculating ellipse. It means that, at a given 
instant of time, the particle may be assumed to lie on the Keplerian ellipse determined by 
the position and velocity at that instant thought as initial conditions for an unperturbed 
motion; at the next instant of time the particle will be at a point of another Keplerian 
ellipse, slightly different with respect to the previous one and determined by the real 
position and velocity of the test body at the new instant of time thought as new initial 
conditions for an unperturbed Keplerian orbit. In other words, the real trajectory of 
the test body at every instant may be approximated by an osculating Keplerian ellipse. 
So the perturbed motion can be described in terms of unperturbed Keplerian elements 
varying in time. See fig. 1. 

Consider the frame K{x, y, z} previously defined and a frame K {X, Y, Z} with 
the Z axis directed along the orbital angular momentum I of the test body, the plane 
{X, Y} coinciding with the orbital plane of the test particle and the X axis directed 
toward the pericenter. K and K have the same origin located in the center of mass of 
the central body. The Keplerian orbital elements are [1] 

• a, e, the semimajor axis and the ellipticity which define the size and the shape of 
the orbit in its plane 
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Fig. 1. - Orbital elements for a general Keplerian osculating ellipse. 



• r2, 2, the longitude of the ascending node and the inchnation which fix the ori- 
entation of the orbit in the space, i.e. of K with respect to K. The longitude of the 
ascending node fl is the angle in the equatorial plane of K{x, y, z} between the x axis 
and the hne of nodes in which the orbital plane intersects the equatorial plane. The 
inclination i is the angle between the z and Z axis 

• CO, Ai, the argument of pericenter and the mean anomaly. The argument of 
pericenter uj is the angle in the orbital plane between the line of nodes and the X axis; it 
defines the orientation of the orbit in its plane. The mean anomaly A4 specifies the motion 
of the test particle on its orbit. It is related to the mean motion n = (GM)^/^a~^/^, 
where M is the mass of the gravitating central source, through A4 = n{t — tp) in which 
tp is the time of pericenter passage 

It is customary to define also 

• w = Q + UJ, the longitude of pericenter 

• u = CO + f, the argument of latitude where / is the angle, called true anomaly, 
which in the orbital plane determines the position of the test particle with respect to the 
pericenter 

• s = w + n{tf) — tp), the mean longitude at the epoch to- If to = 0, it is customary 
to write e as Lq — w — ntp. 

The matrix R^x for the change of coordinates from K to ii' is [1] 

(cosficosw — sin f2 COS i sin (J — cosfisinw — sin f2 cos « cos w 
sin O cos CO + cos Cl cos i sin cj — sin f2 sin a; + cos Cl cos i cos co 
sin i sin co sin i cos u) 

Using eq. (16) and X = rcos/, Y = rsmf, Z = it is possible to express the 
geocentric rectangular Cartesian coordinates of the orbiter in terms of its Keplerian 



sin n sin i 
- cos O sin i 
cosi 
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elements 

X = r(cos u cos r2 — sin u cos i sin fi) 

(17) ^ y = r(cos wsinfi + sinwcoszcosfi) 

z = rsinusini. 

Redefining suitably the origin of the angle O so that cosfi = 1, sinO = 0, the previous 
equations become 

: r cos u 

(18) ^ y = rsinwcos?' 

: rsin?isin2. 

Considering for the test particle the total mechanical energy with the sign reversed, 
according to [27], = —Stot — ^{T + ^)j where T and hi are the kinetic and potential 
energies per unit mass, it is possible to work out the analytical expressions for the rate of 
changes of a, e, i, fl, u, M. due to any non-central gravitational contribution. To this 
aim it is useful to isolate in U the central part —C of the gravitational field from the terms 
—TZ which may cause the Keplerian orbital elements to change in time: U = —C — TZ. 
For a spherically symmetric body, becomes 

^ GM ^ ^ GM ^ 
r 2a 

Concerning the pcrturbativc scheme to be employed, the well known Lagrange planetary 
equations are adopted. At first order, they are 



(20) 
(21) 



da 2 dTZ 



dt na dM 

de l-e^ dTZ (1-6^)1/2 gji 



dt na^e dM. na'^e duj ' 

dd _ . 1 dTZ 1 dTZ 

dt na2(l — e2)i/2 sini du) na2(l — e^)!/^ sini 50' 

dn 1 dTZ 



(23) 



dt na2(l — 6^)1/2 sin i di ' 



div . 1 dTZ (l-e2)V2 dTZ 

= ~<=°^*^71 2M/2 • ■ 1^ + 2 -TT' 

dt na2(l — e2)^/2 sm^ di na'^e de 

dM l-e"^ dTZ 2 dTZ 



(25) 



dt na^e de na da 



The idea consists of using Cg^ to obtain a suitable non-central term TZgm to be employed 
in these equations. This can be done considering the Hamiltonian for the test particle 

(26) W = p-v-£. 
Inserting eq. (13) in eq. (26) one has 

(27) W = WW+Wg„, 
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with Hgm = —^{-^g • v). So it can be posed 

(28) T^gm 



Hgrn __ (Ag ■ V) 



Now it is useful to express eq. (28) in terms of the Keplerian elements. Referring to 
cq. (9), eq. (18), and recalling that in the frame K{x, y, z} J — (0, 0, J) and that for 
an unperturbed Keplerian motion 

(29) = 

r a{l — e^) 

it is possible to obtain, for a perfectly spherical central body 



2GJcosi. 2GJcosi (1 + ecos/) . 

a(l-e2) 



(30) = --^—-U = o JJ u. 



In eq. (30) u ~ / is assumed due to the fact that the osculating element a; may be 
retained almost constant on the temporal scale of variation of the true anomaly of the 
test body. 

2'4. Secular gravitomagnetic effects on the Keplerian orbital elements. - The secular 
effects can be worked out by adopting the same strategy followed, e.g., in [27] for a similar 
kind of perturbing functions. By using averaging we implicitly make several assumptions 
about the system that may restrict the validity and applicability of our result. Averaging 
is generally applied to systems where the perturbing force is sufficiently small so that, over 
one orbital period, the deviations of the true trajectory from the Keplerian trajectory 
are relatively small. This perturbation can generally be related to the magnitude of 
the perturbing acceleration as compared to the central-body attraction. In this case, 
clearly, wc can view the application of the averaging approach as valid. When cq. (30) is 
mediated over one orbital period of the test body, a, e, i, f2 and u) are to be considered 
constant 

~ G2J _ ^G2Jcosi r^^df __^G 2Jcosi 



1 r"^ G2J . , G2Jcosi f 

(31) (^)2.=-py^ ^V'°'"^''=-''^^^i 



a(l- 



The relation du = du + df = df has been used in eq. (31). Eq. (31) can now be used in 
determining the secular changes of the Keplerian orbital elements of the test body. Prom 
it one gets 

,33, ^ = 0, 

(34, 2<^=0, 



(35) 
(36) 



9(7?.) 2^ G 2nJsmi 
di (P- a(l — e^) ' 

G 4nJecosi 
de c^a(l — 6^)2' 
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A particular care is needed for the treatment of n when the derivative of (7?-) 2^ with 
respect to a is taken; indeed, it must be posed as 

dn 



In eq. (37) 
and 

dn G 3t?, Jcosz 
^ da (? a2(l — e^) ' 

From cq. (20) and eq. (32) it appears that there are no scciilar changes in the semimajor 
axis(^), and so the orbital period of the test body, related to the mean motion by P = 
27r/n, can bo considered constant. So in eq. (37) only eq. (38) must be retained. Using 
eqs. (32)- (37) in eqs. (20)- (25) one obtains for the averaged rates 

0, 
0, 

0, 

G 2J 
?a3(l-e2)3/2' 

G 6 J cos i 
~?a3(l-e2)3/2' 

0. 

They arc the well known Lcnsc-Thirring equations [7, 8, 25]. In deriving them it has 
been assumed that the spatial average over / yields the same results for the time average 
[33, 34, 2]. 

The effects of the quadrupole mass moments of the source arc treated in [25]: it turns 
out that they are negligible in the context of the LAGEOS Lense-Thirring experiment. 

3. — Orbital perturbations induced by solid Earth and ocecin tides 

3'1. Introduction. - In this section the problem of calculating the gravitational time- 
dependent perturbations induced by solid Earth and ocean tides on the Keplerian orbital 



(*) Here only the effect of the off-diagonal gravitomagnetic components of the metric tensor is 
considered, neglecting any other physical influence of gravitational and non gravitational origin. 



(37) 



d{n),^ d{n) 



2-K 



da 



da 



+ 



dn 



2-K 



G 2nJcosi 
(? o?{\ — e^) ' 



(39) 



2tt 



dn 



(40) 
(41) 

(42) 
(43) 
(44) 

(45) 



da 
H 
de 
di 
di 
di 

dn 

Itt 
duj 
'dt 
dM 



LT 



LT 



LT 



LT 



LT 



LT 
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elements of LAGEOS and LAGEOS II is addressed. The focus of the analysis is on the 
nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II in view of their close 
connection to the measurement of the Lcnsc-Thirring effect. The pcrturbative scheme 
adopted is that based on the Lagrange' s planetary equations [1, 35, 36, 37]. 

An evaluation from first principles of the amplitudes, the periods and the initial phases 
of the tidal perturbations on the nodes of LAGEOS and LAGEOS II and the perigee of 
LAGEOS II is useful, in the present context, for the following reasons. 

• It allows one to point out which tidal constituents the Lense-Thirring shift is mainly 
sensitive to, so that people can focus the researches on them 

• It allows one to get insight in how to update the orbit determination softwares. In 
this way the impact on the time series of those tidal constituents which will turn out 
to be more effective could be reduced along with the number of the parameters to be 
included in the least-squares solution 

• For a given observational time span Tobs it allows one to find those constituents 
whose periods are longer than it and consequently may act like superimposed bias (e.g. 
the 18.6-year and 9.3-year zonal tides and the Ki I = 3 p = 1 q = —1 ocean diurnal 
tide on the perigee of LAGEOS II) corrupting the detection of the gravitomagnotic shift. 
In these cases, if we know their amplitudes, periods and initial phases we could fit and 
remove them from the time series(^) or, at least, it should be possible to assess the 
mismodelling level induced by them on the detection of the Lense-Thirring trend. So, 
also these estimates can be included in the final error budget. (See section 4) 

• It can be used as starting point for numerical simulations of the combined residual 
data in order to check, e. g., the impact of the diurnal (to=1) and semidiurnal (m=2) 
tidal perturbations (not cancelled out by the Ciufolini' s observable), as done in section 
4 

The section is organized as follows. In subsect. 2 the essential features of the tide 
generating potential, from which the tidal forces are traditionally derived, are reviewed. 

The solid Earth tidal perturbations are worked out in subsect. 3. An evaluation 
of their mismodelling and a comparison with the Lense-Thirring perturbations on the 
nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II over 4 years are also 
presented. 

Subsect. 4-5 are devoted to the orbital perturbations due to the ocean tides and their 
mismodelling in connection with the Lense-Thirring drag. They are treated in detail 
because the ocean tides are known less accurately than the solid ones, so that the impact 
of their mismodelling on the gravitomagnetic trend is more relevant. 

Subsect. 6 is devoted to the conclusions. 

3"2. The tide generating potential. - On the scale of the Earth's dimensions the grav- 
itational field of any not too far astronomical body B cannot be considered uniform so 
that the various points of our planet are acted upon by a differential gravitational pull 
in the external field of B. This is the origin of the Earth's solid, ocean and atmospheric 
tides, mainly duo to the Moon and the Sun. It is customary to derive the tidal forces 
acting on a point on the Earth' s surface from the so called tide generating potential 



(^) If the period of the disturbing harmonic, assumed to be known, is shorter than the time 
scries length the perturbation can be viewed as an empirically fit quantity. But if its period is 
longer than the time series length it is not possible to fit and remove the harmonic from the 
data without removing also the true trend, unless we know the initial phase. 
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which, in its most general form, for a perturbing body B can be written as [38, 39, 20] 

GMb ^ f r y 



3 R"^ / r 



3 r 



4 4 r 

-P2(cos2;) + 3^-P3(cos2;) 



In cq. (46) the following quantities appear 

• R mean Earth' s equatorial radius 

• Mb mass of the perturbing body 

• (Ib instantaneous distance between the Earth and the perturbing body B 

• r distance between the center of mass of the Earth and a point on its surface 

• Pi{x) Legendre polynomials of degree I 



(47) 



(-1)' s 

l\2' 



(1 - 



• z geocentric zenithal distance of the perturbing body 

• Cb mean distance between the Earth and the perturbing body 
The factor 



(48) 

is often rewritten as 



-GMb^(^) 
4 c% \rJ 



(49) 
with 



DBir) 



Mb 



Mm c 



D{r), 



(50) 



in which 



Dir) 



3 „ , r R'^ 

-GMm^ 



(51) 



3 7?^ 
£)i = - GMm — 

4 C^j 



is the Doodson constant for the Moon, having the dimension of an energy per imit mass, 
and Mm is the Moon's mass. Thus the tide generating potential outside of the Earth $^ 
is a harmonic function of degree I and can be written in the form $^ = r''Yi{(f>, A) where 
Yi((t), A) is a surface harmonic, with (j) and A the terrestrial latitude and East longitude. 
At the Earth's surface it becomes $^(1^) = 

It is useful to express eq. (46) in terms of the geocentric equatorial coordinates 
{(/), 6b : Hb} of an Earth- fixed frame where 

• Sb declination of the perturbing body 

• Hb hour angle of the perturbing body 
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This can be done by means of the foUowing spherical trigonometric formula 

(52) cos^; = sin sin + cos (f) cos 5b cos Hb- 

The result is the Laplace's tidal development: it consists of a sum of terms each of which 
is the product of a factor depending only on the latitude of a given place on the Earth 
(the geodetic coefRcient) and a time-dependent factor which depends on the astronomical 
coordinates 6b and Hb of the perturbing body. E.g., the term of degree I = 2 becomes 

[(i^)(i-.„^..). 

(53) + sin 2(j) sin 25 b cos Hb + cos^ (f) cos^ 5 b cos 2Hb ] • 

In this sum the quantities CB/ds, 5b, Hb exhibit a very complex behavior in time due 
to the astronomical motions of the perturbing body B. They must be carefully expanded 
in periodic components in order to obtain an expression for the tidal potential as a sum 
of harmonic terms. This can be done by using the ephemerides tables for the perturbing 
bodies. Obviously, the most relevant arc the Moon and the Sun. Since the ephemerides 
refer to the ecliptic and not to the celestial equator, it is necessary to express 5b and Hb 
in terms of ecliptical coordinates. Spherical trigonometry gives 

(54) sin (5b = sinecos/Ss sin As + cosesinAs, 

(55) cos5bcosHb = cos (3 b cos Xb cos 9 + (cose cos /3b sin A^ — sin e sin /S^) sin 0, 
where 

• e inclination of the ecliptic to the celestial equator 

• As , Pb ecliptical longitude and latitude of the perturbing body, respectively 

• 9 sidereal time 

The lunar and solar ephemerides express rfs. As, ,9s in terms of time-dependent 
series of sines and cosines whose argiimcnts arc linear combinations of the mean Lunar 
time T and other astronomical quantities of the Moon and the Sun (®). Among them the 
most useful are the mean ecliptical lunisolar variables s, h, p, N , [38]. They are 

• s mean ecliptical longitude of the Moon with period P, = 27.32 days 

• h mean ecliptical longitude of the Sun with period = 365.2422 days 

• p mean ecliptical longitude of the Moon's perigee with period Pp = 3, 232 days 
(8.84 years) 

• N = —N mean ecliptical longitude of the Moon's node with the sign reversed, 
with period Pn = 6,798.38 days (18.61 years). The lunar node, moving in the Earth- 
fixed frame from East to West counterclockwise, is retrograde with respect to the other 
lunisolar longitudes 

• Ps mean ecliptical longitude of the perihelion with period Pp^ = 7.65 x 10^ days 
(20,953 years) 



(^) The most recent developments take in account also the influence of other planets like Venus 
and Jupiter [47]. 
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Putting eqs. (54)-(55) in the Laplace' s tidal development, with the ephemerides 
expansions for rfs, Ab, (3b, allows one to obtain the harmonic expansion of the tide 
generating potential $b [40, 41, 42, 43, 44, 45, 46, 47]. It consists of a sum of terms, the 
constituents, with sinusoidal time dependence where the sines and cosines have arguments 
involving linear combinations at of the previously defined orbital longitudes of the Sun 
and the Moon 

(56) at = jiT + j2S + j^h + jip + jr,N + jePs- 

The circular frequency of the tidal bulge generated by the corresponding constituent, 
viewed in the terrestrial frame, is then given by 

(57) a = 211 f = jiT + j2S + j-sh + j^p + j^N + jePs] 
by considering that t = — s, eq. (57) becomes 

(58) a = ji6 + (j2 - ji)s + jsh + jip + j^N + j^p^. 

The coefficients jk, fc = 1,...,6 are small integers which can assume negative, positive 
or null values. The advantage of using such lunisolar ecliptical variables relics in the 
fact that they, over an interval of a century or so, are practically linear increasing with 
time. This feature will reveal itself very useful in integrating the equations for the 
orbital perturbations. Each tidal constituent is identified by the set of the six integers 
jk, /c = 1, 6 arranged in the so called Doodson number 

(59) ii {32 + 5) (is + 5) . (i4 + 5) (is + 5) (je + 5) . 

In it the integer j\ classifies the tides in long period or zonal {j\ = 0), diurnal or tesseral 
[ji = 1) and semidiurnal or sectorial (ji = 2). 

The tide generating potential is the cornerstone for every calculation of tidal per- 
turbations on the satellites' Keplerian orbital elements. Indeed, if the Earth, globally 
modeled as a nonrigid body, is acted upon by a tidal potential harmonic in degree, it 
deforms itself giving rise, among other things, to a periodic redistribution of the masses 
within its volume. This deformation acts upon a point in the external space surrounding 
the Earth by means of an additional potential [48] 

+ l /n^2^+l 

(60) 0;(r) = fcjfij ^i{H) = ki(j] $;(r) 

where ki is one of the so called Love numbers (they will be defined precisely in the 

following sections) . A close orbiting satellite is perturbed by such an additional potential 
and senses the global effect of the solid and fluid mass redistribution in the Earth. 
This means that if one wish to employ a pool of satellites in order to recover the tidal 
parameters, like the Love numbers, entering the tidal perturbations and wants to use the 
so obtained values to predict the perturbations on some particular satellite, he or she 
must be very careful. Indeed, such values are aliased by the whole of the effects sensed 
by the measuring satellites; it is as if they would be "tailored" for the satellites used in 
their recovery. So, these "effective" tidal parameters neither can be directly compared 
to those measured on the Earth's surface or predicted by any theoretical model and 
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related to the physical properties of the Earth, nor can be used in calculating a priori 
the tidal perturbations on the satellites of interest for some particular application, as 
is the case for the detection of the Lonsc-Thirring effect. Differently stated, when the 
tidal perturbations are to be predicted for a given satellite, it is incorrect to use for it 
the "effective" values of the geophysical parameters recovered by other satellites. It is so 
customary to split the analytical expressions of the tidal perturbations into a part due 
to the solid Earth tides and another one which accounts for the effects of the oceans [49] . 
After this step, their effects on the satellite of interest must be analytically predicted 
using for the various parameters the values theoretically calculated or measured in such 
a way that they reflect the effective Earth's properties. 

3'3. Solid Earth tidal perturbations. - Concerning the solid Earth tides [50], as in 
[35, 51], the starting point is the frequency-dependent model of [52]. It is based on the 
assumption of a perfectly elastic, hydrostatically prestressed, ellipsoidal rotating Earth 
acted upon by the lunisolar tidal potential. The interior of our planet is thought to 
be made of a solid inner core, a fluid outer core, and a solid mantle capped by a thin 
continental crust, without oceans and atmosphere. Substantially, the Earth is thought 
as a set of coupled harmonic oscillators showing a variety of resonant frequencies, the 
normal modes, which can be excited by the external forcing constituents of the tide 
generating potential [53]. 

The expansion adopted for it is that of Cartwrigt and Tayler in [41] in terms of 
harmonic constituents of degree I, order m and circular frequency a 



Only the real part of eq. (61) has to be retained. The quantities entering eq. (61) are 
• Y™ spherical harmonics 



in which the Condon-Shortley phase factor (—1)"* has been neglected in order to compare 
our results to those of [51] 

• Hj^ (/) coefficients of the harmonic expansion. They contain the lunisolar ephemerides 

information and define the modulus of the vertical shift V/ g in the equipotential surface 
at the Earth's surface for r = with respect to the case in which only proper Earth's 
gravity is considered. The values used in the present calculation are those recently re- 
leased in [47] using the ELP2000-85 series for the ecliptical coordinates of the Moon 
[54] and the VSOP87 series for those of the Sun [55]. Their accuracy is of the order of 
10~^ m. In order to use them in place of the coefficients of [42], which have a different 
normalization, a multiplication for a suitable conversion factor [56] has been performed. 
Due to the extreme smallness of the ff["(/) for / > 3, in the following summations only 
the term of degree I = 2 are to be considered 

• Cim additive constant. It is equal to — 7r/2 in order to generate sines for the 
I — m odd constituents, while it is equal to in order to yield cosines for the I — m even 
constituents. 



(61) 




(62) 
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The Earth free space potential under the action of the tidal constituent of degree I, 
order m and circular frequency a, dropping the time dependence due to e"^*, can be 
written as 



(63) 



= Re gH[' 



JO) 



R 



i+i 



-I yr+K^[- 



R\ 



i+3 



W+2 



in which 

• r distance from the Earth' s center of mass and a point in the space outside of the 
Earth 

• The parameter 



(64) 



is the modulus of the Love number k [48, 57] for a static, spherical Earth. In general, it is 
defined as the ratio of the free space gravitational potential (/'(r), evaluated at the Earth' s 
equator, and the tide generating potential $(r) calculated at the mean equatorial radius. 
If the rotation and non-sphericity of the Earth are accounted for, there are different 
values of k for any l,m and /. In general, they are evaluated by convolving, in the 
frequency domain, the tide generating potential with the transfer function for a non 
rigid Earth [58, 59]. There are several theoretical calculations for the Love numbers k 
[52, 60, 57, 56, 61]. They differ in the choice of the Earth's interior models [77, 63] 
adopted for the calculation of the transfer function and the departure from elasticity of 
the mantle' s behavior. Particularly important is also if they account for, or not, the 
normal modes in the diurnal band. In eq. (64) the lERS conventions [56] have been used 



(65) 



kUf) 



Imkim + Skf^^, 



(66) kt,{f) = Rekim + 5kf, 

where Im kim and Re kim are the frequency-independent parts of Love numbers, and 
Skj"^'^'^, Skj- are the frequency-dependent corrections. They are important in the diurnal 

band, through Skj^, and in the zonal band with Jfc^"*^' due to the anelasticity of the 
mantle. The k'l^{f) account for the rotation and ellipticity of our planet. In the present 
analysis for the latitude dependent Love numbers k^{f ) the values quoted in [61] have 
been adopted. 

Neglecting, in this first step, their tiny contribution, the free space potential by means 
of which the Earth responds to the forcing lunisolar tidal potential can be written, in 
the time domain, as 

Here and in the following, the upper expressions refer to l — m even while the lower ones 



cos {(Tt + mX) 
sin {at + m\) 
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refer to I — m odd. Eq. (67) can be put into the form 



(68) 



GM ^ 1 /R^ ' 



EE 



1=2 m=0 



{Cim COS mX + Sim sinmA)P™(sin 0), 



in which, assuming temporarily for the sake of simphcity that the Love numbers are not 
frequency-dependent, the coefficients Cim and Sim are 



(69) 
(70) 



Cim = Aimkl'2^Hl 



Sim — ^/m^/TO y^-^r 
/ 



cos at 
sin at 

— sin at 
cos at 



with 



(71) 



Air, 



l2l + l{l-m)\ 
4n (Z + rn)!' 



The Cim and Sim have the dimensions of lengths. Eq. (68) is formally equal to the well 
known expression of the static geopotential worked out by Kaula in ref.[l] 



(72) 



GM 



1=2 m=0 

By means of the transformation 



1 + ^ ^ f — j Pim (sin (j)) [Cim COS mX + Sim sin mX] 



(73) i Prism<f>)e'^^ = i f^Xj^PimA^) E ^'^^(e) 

^ ^ ^ ^ p—O q— — oo 



with 

(74) i^impq = (/ - 2p)uj + {l-2p + q)M + m{D. 

it can be cast into the familiar form 



(75) 



GM - , ^ ]^ f R ^ ^ 

~R~EEa("^) E E -^impl*) Gipg(e)5',™pg(w, ri,7W, 0), 



1=2 m=0 



p=0 q=—oo 



where Fimp{i) and Gipq{e) are the inclination and eccentricity functions respectively [1]. 
The Simpq is given by 



(76) 



Si 



mpq 



Cim 

—Si„ 



COS i^lmpq + 



Si„ 

Clr. 



sin tpi 



mpq- 
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Eq. (76) in our case becomes 
(77) 



Slmpq = ^lmk[^ ^ H]^ COS {at + Iplmpq)- 
f 



It is worthwhile noting that eq. (73) expresses a transformation of coordinates from 
the geocentric equatorial frame rotating with the Earth to the geocentric inertial frame 
K{x, y, z}. This fact will reveal itself of paramount importance in evaluating the 
frequencies of the tidal perturbations on the Keplerian elements. 

Up to now the response of the Earth to the forcing tidal potential has been considered 
as perfectly elastic; if it had been so, there would be no delay between the Earth's tidal 
bulge and the Moon in the sense that when the latter passes at the observer's meridian, 
say, at A the tidal bulge reaches A at exactly the same time. The reality is quite 
different since complicated mechanisms of energy dissipation in the interior of the Earth 
[64] makes its response to depart from the perfectly elastic behavior previously sketched. 
A phase lag with respect to the external lunisolar potential must be introduced in the 
sense that the tidal bulge reaches the observer's meridian at A after a certain time At 
with respect to the passage at A of Moon. It is generally assumed that it is the Earth's 
mantle to exhibit an anclastic behavior at the various frequencies of the tide generating 
potential; this topic is not yet well understood, but if one wants to account for the phase 
lag introduced by it, he or she has to adopt complex frequency-dependent Love numbers 
[57]. In order to give a quantitative estimation also of the lag angle it has been decided 
to adopt the values of ref.[57] quoted in ref.[56]; indeed, the authors of ref.[61] do not 



quote the imaginary part of (/) and the discrepancies with their values for k^!^ (/) 

amount to only 0.1 %. Moreover, concerning the effects of mantle's anelasticity, the 

results of ref.[61] arc not less uncertain than the other ones. 

An equivalent form of eq. (69) and eq. (70), for a single constituent of frequency a, 
which reveals useful in handling with complex quantities is 



,(0)/ 



(78) 



(/)e^' 



with T] = —i if I — m is odd and r/ = 1 if Z — m is even. If the anelasticity of the Earth's 
mantle is to be accounted for writing kj^{f) = ki^{f) + iki^{f), eq. (78) becomes 

Cimf - iSimf = r]AimHp{f) X [(fc(^ COS cri - k{^ sin at) + 



(79) 



or 



+i{khn sin ai + kf^ cos at)] 



(80) 
(81) 



Clr, 



Sin 



Aim X] 

f 

Aim. > HP 



f 



^i^coso-t-fcf^sinat 
^/mSinCTt + fcf^coscTi 

-khn sin at -k{^ cos at 
^ZmCOSCTt-fcf^sincTi 
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If eq. (80) and eq. (81) are put into eq. (76), it is straightforward to obtain for it 

(82) Slmpg = Ai^ ^r^l™ (/) cos {at + i^lmpg - Slmf), 

f 

with 

(83) if) = ^/kKiff + kUff, 



(84) tan^im/ = 



kLifY 



The factor Simf is the phase lag of the response of the sohd Earth with respect to the 
tidal constituent of degree I, order m and circular frequency a induced by the anelasticity 
in the mantle: notice that if ki^{f) = 0, also ta,n Simf = 0. An inspection of the values 
quoted in [56] shows that, at the low frequencies of the zonal constituents of the tide 
generating potential, the role played by the mantle's anelasticity is more relevant than 
in the diurnal and semidiurnal bands. Indeed, while for the tesseral and sectorial tides 
one finds = —0.00144 and fc^j = —0.00130 for all the frequencies, in the zonal band 
Hoif) varies from -0.00541 to -0.00192. 

Eq. (75), with the Simpq given by eq. (82), is the dynamical, non-central part of the 
geopotential due to the response of the solid Earth to the forcing hmisolar tidal action. 
In the linear Lagrange equations of perturbation theory for the rates of change of the 
Keplerian elements it can play the role of the perturbative term TZ. It is an observed 
fact that the secular motions are the dominant perturbation in the elements oj, J7, M. 
of geodetically useful satellites. So, taking as constants a, e, i and consider as linearly 
variable in time uj, fl, A4 and 6, apart from at, the following expressions, at first order, 
may be worked out 



Anf = — ^J= — yy - 



I +00 ^ 



(85) 



p=0 g=— oo 



Au)f 



EE(f 

1=0 m=0 ^ 



/ +00 r 



(86) xAr^Y E 

p=0 g= — oo 

with 



1 — „ dGipq cos i dFi 



Imp 



Imp 



de sin i di 



G 



Ipq 



(87) jfimpq = {l- 2p)a; + {l-2p + q)M + m(Q - 6) + at - Si^f, 
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(88) fp = {l- 2p)w + {l-2p + q)M + m{tl -e) + a. 

The indirect, second order effects due to the oblateness of the Earth [65] have been 
neglected in this section. Due to the secular trend of the Lense-Thirring effect, only the 
perturbations whose periods are much longer than those of the orbital satellites motions, 
which, typically, amount to a few hours, are to be considered. This implies that in such 
terms the rate of the mean anomaly does not appear and the condition 

(89) l-2p + q = Q 

must hold. Moreover, if the effect of Earth's diurnal rotation, which could introduce 
periodicities of the order of 24 hours, is to be neglected, one must retain only those 
terms in which the non-negative multiplicative coefficient ji of the Greenwich sidereal 
time in at coincidc's to the order m of the tidal constituent considered: in this way in 
fp the contributions of 9 are equal and opposite, and cancel out. With these bounds on 
Z, m, p and q the circular frequencies of the perturbations of interest become 

(90) fp = tf + {l-2p)L0 + mh 
in which 

(91) tf = {j2 - m)s + js A + jip + j^N + jeps . 

In the performed calculation only the degree I = 2 constituents have been considered due 
to the smallncss of the and H^[f). For I ~ 2, p runs from to 2, and so, in virtue 
of the condition I — 2p + q = 0, q assumes the values —2, 0, 2. From an inspection of 
the table of the eccentricity functions Gipq{e) in [1] it turns out that G20-2 = G222 = 0, 
while G210 = (1 — e^)~^/^. For this combination of I, p and q the condition I — 2p = Q 
holds: the frequencies of the perturbations are, in this case, given by 

(92) fp = tf+mQ. 

Passing from a to Tf + mil, as previously noticed, is equivalent to leaving the Earth- 
fixed frame of reference for the inertial one, in which the relevant physical feature is 
the orientation of the tidal bulges with respect to the orbital plane of satellite. Indeed, 
T f+mCl could be considered as the relative frequency of motions of the tidal bulge and the 
orbital plane, both viewed in the inertial frame. T/ is the opposite of the inertial rate of 
the tidal bulge obtained from the Earth-fixed one subtracting the Earth's rotation rate 9 
times m in order to account for the symmetry of the bulges. The periodicities introduced 
in the perturbations by eq. (92) are, in general, much longer that the ones present in 
the tidal potential because, while the latter are dominated by ji9 due to the rotation of 
the Earth-fixed frame, the former are mainly determined by the precessional rates mil 
of the orbital frames of the different satellites in the inertial frame. This explain the 
major variety of periodicities in the orbital elements' perturbations of the near-Earth's 
satellites with respect to the tides sensed by an observer on the Earth surface which 
are mainly concentrated in the diurnal and semidiurnal bands. In eqs. (85)-(86) it is 
worthwhile noting that the amplitudes of the perturbations are inversely proportional to 
the frequency of the perturbation fp. This implies that some tidal constituent which in 
the Earth-fixed frame has an high frequency, in the inertial frame may induce relevant 
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perturbations because its frequency greatly decreases; this fact could compensate an 
eventual small value of its coefficient in the harmonic expansion of the tide generating 
potential (calculated in the terrestrial frame). In other words, constituents which on 
the Earth's surface would produce small tides may become important in perturbing the 
orbits of near-Earth' s satellites due to the changes in their frequencies when viewed in 
the inertial frame. 

3'4. Discussion of the numerical results. In view of their direct implication in eq. 
(5), in the following we shall focus our attention to the nodes of LAGEOS and LAGEOS 
II and the perigee of LAGEOS II. In Table I, Table II, and Table III the results for them 
are shown; since the observable quantity for cj is ead) [14], the calculation for the perigee 
of LAGEOS, due to the notable smallness of the eccentricity of its orbit, have not been 
performed. 

Table I. - Perturbative amplitudes on the node f2 of LAGEOS due to solid Earth tides for 
I = 2, p = 1, q = 0. In the first column the Doodson number of each constituent is quoted 
followed by the Darwin's name, when it is present. The tidal lines are listed in order of decreasing 
periods. The coefficients HJ^{f) are those recently calculated by Roosbeek and multiplied by 
suitable normalization factors (lERS standards) in order to make possible a comparison with 
those of Cartwright and Edden. tan Simf expresses the phase lag of the solid Earth response with 
respect to the tidal potential due to the anelasticity in the mantle. The periods are in days, the 
amplitudes in mas and the Hp{f) in meters. 



Tide 




Period 


Amplitude 


fca Love number 


Hrif) 


tan Simf 


055.565 




6,798.38 


-1,079.38 


0.315 


0.02792 


-0.01715 


055.575 




3,399.19 


5.23 


0.313 


0.000272 


-0.015584 


056.554 


Sa 


365.27 


9.96 


0.307 


-0.00492 


-0.01135 


057.555 


Ssa 


182.62 


31.21 


0.305 


-0.03099 


-0.01029 


065.455 


Mm 


27.55 


5.28 


0.302 


-0.03518 


-0.00782 


075.555 


Mf 


13.66 


4.94 


0.301 


-0.06659 


-0.007059 


165.545 




1232.94 


-41.15 


0.259 


-0.007295 


-0.00554 


165.555 


Ki 


1,043.67 


1,744.38 


0.257 


0.3687012 


-0.0055933 


165.565 




904.77 


203.02 


0.254 


0.050028 


-0.005653 


163.555 


Pi 


-221.35 


136.44 


0.286 


-0.12198 


-0.005017 


145.555 


Oi 


-13.84 


19 


0.297 


-0.26214 


-0.00484 


135.655 


Qi 


-9.21 


2.42 


0.297 


-0.05019 


-0.00483 


274.556 




-1,217.55 


1.68 


0.301 


0.000625 


-0.00431 


274.554 




-1,216.73 


-6.63 


0.301 


-0.00246 


-0.00431 


275.555 


K2 


521.835 


-92.37 


0.301 


0.0799155 


-0.00431 


273.555 


S2 


-280.93 


182.96 


0.301 


0.2940 


-0.00431 


272.556 


Ta 


-158.80 


6.04 


0.301 


0.0171884 


-0.00431 


255.555 


Ma 


-14.02 


19.63 


0.301 


0.6319 


-0.00431 


245.655 




-9.29 


2.49 


0.301 


0.12099 


-0.00431 
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Table 11. - Perturbative amplitudes on the node of LAGEOS II due to solid Earth tides 
for I = 2, p — 1, q = 0. In the first column the Doodson number of each constituent is quoted 
followed by the Darwin's name, when it is present. The tidal tines are listed in order of decreasing 
periods. The coefficients H]^{f) are those recently calculated by Roosbeek and multiplied by 
suitable normalization factors (lERS standards) in order to make possible a comparison with 
those of Cartwright and Edden. tan Simf expresses the phase lag of the solid Earth response with 
respect to the tidal potential due to the anelasttcity m the mantle. The periods are in days, the 
amplitudes in mas and the Hp{f) in meters. 



Tide 




Period 


Amplitude 


k2 Love number 


Hr{f) 


tan Simf 


055.565 




6,798.38 


1,982.16 


0.315 


0.02792 


-0.01715 


055.575 




3,399.19 


-9.61 


0.313 


-0.000272 


-0.015584 


056.554 


Sa 


365.27 


-18.28 


0.307 


-0.00492 


-0.01135 


057.555 




182.62 


-57.31 


0.305 


-0.03099 


-0.01029 


065.455 


Mm 


27.55 


-9.71 


0.302 


-0.03518 


-0.00782 


075.555 


Mf 


13.66 


-9,08 


0.301 


-0.06659 


-0.007059 


165.565 




-621.22 


-58.31 


0.254 


0.050028 


-0.005653 


165.555 


Ki 


-569.21 


-398 


0.257 


0.3687012 


-0.0055933 


165.545 




-525.23 


7.33 


0.259 


-0.007295 


-0.005541 


163.555 


Pi 


-138.26 


35.65 


0.286 


-0.1219 


-0.005017 


145.555 


Oi 


-13.33 


7.66 


0.297 


-0.26214 


-0.00484 


135.655 


Qi 


-8.98 


0.98 


0.297 


-0.05019 


-0.00483 


275.555 


K2 


-284.6 


-92.51 


0.301 


0.079915 


-0.004318 


274.556 




-159.96 


-0.40 


0.301 


0.000625 


-0.00431 


274.554 




-159.95 


1.6 


0.301 


-0.00246 


-0.0043 


273.555 


52 


-111.24 


-133.04 


0.301 


0.29402 


-0.00431 


272.556 


T2 


-85.27 


-5.96 


0.301 


0.017188 


-0.00431 


255.555 


M2 


-13.03 


-33.05 


0.301 


0.6319 


-0.004318 


245.655 


N2 


-8.84 


-4.35 


0.301 


0.12099 


-0.00431 
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Table III. - Periurbative amplitudes on the perigee to of LAGEOS II due to solid Earth tides 
for 1 = 2, p = 1, q = 0. In the first column the Doodson number of each constituent is quoted 
followed by the Darwin's name, when it is present. The tidal lines are listed in order of decreasing 
periods. The coefficients HP{f) are those recently calculated by Roosbeek and multiplied by 
suitable normalization factors (lERS standards) in order to make possible a covipanson with 
those of Cartwright and Edden. tan Simf expresses the phase lag of the solid Earth response with 
respect to the tidal potential due to the anelasticity in the mantle. The periods are in days, the 
amplitudes in mas and the Hp{f) in meters. 



Tide 




Period 


Amplitude 


k2 Love number 


Wif) 


tan 5ijnf 


055.565 




6,798.38 


-1,375.58 


0.315 


0.02792 


-0.01715 


055.575 




3,399.19 


6.66 


0.313 


-0.000272 


-0.015584 


056.554 


Sa 


365.27 


12.69 


0.307 


-0.00492 


-0.01135 


057.555 


Ssa 


182.62 


39.77 


0.305 


-0.03099 


-0.01029 


065.455 


Mm 


27.55 


6.74 


0.302 


-0.03518 


-0.00782 


075.555 


Mf 


13.66 


6.30 


0.301 


-0.06659 


-0.007059 


165.565 




-621.22 


290.43 


0.254 


0.050028 


-0.005653 


165.555 


Ki 


-569.21 


1,982.14 


0.257 


0.3687012 


-0.0055933 


165.545 




-525.23 


-36.52 


0.259 


-0.007295 


-0.005541 


163.555 


Pi 


-138.26 


-177.56 


0.286 


-0.1219 


-0.005017 


145.555 


Oi 


-13.33 


-38.16 


0.297 


-0.26214 


-0.00484 


135.655 


Qi 


-8.98 


-4.92 


0.297 


-0.05019 


-0.00483 


275.555 




-284.6 


-88.19 


0.301 


0.079915 


-0.004318 


274.556 




-159.96 


-0.38 


0.301 


0.000625 


-0.00431 


274.554 




-159.95 


1.52 


0.301 


-0.00246 


-0.0043 


273.555 


S2 


-111.24 


-126.83 


0.301 


0.29402 


-0.00431 


272.556 


Ta 


-85.27 


-5.68 


0.301 


0.017188 


-0.00431 


255.555 


Ma 


-13.03 


-31.9 


0.301 


0.6319 


-0.004318 


245.655 


N2 


-8.84 


-4.15 


0.301 


0.12099 


-0.00431 



The tidal lines for which the analysis was performed have been chosen also in order 

to make a comparison with the pcrturbativc amplitudes of the ocean tides based on the 
results of EGM96 gravity model [19] which will be shown in the next sections. It is 
interesting to notice that the periods of many of these tidal perturbations are almost 
equal to or longer than an year, a time interval in which their effect may alias the Lense- 
Thirring effect, if the observations are taken on a such temporal scale; indeed, only if the 
data were acquired on time intervals of the order of many years it should be possible to 
clearly resolve the action of these tidal lines with respect to the secular Lcnsc-Thirring 
trend. And also in this case, it could remain the aliasing effect of the 18.6-year and 
9.3-year tides. Up to now, in refs.[66, 15] the authors have analyzed the data acquired on 
a 4 years interval for a suitable combination of residuals showing that many tidal signals 
with periods shorter than it can be resolved by means of suitable fits. 

The results presented in Table I, Table II and Table III can be compared to those of 
[35, 51]. See also [67]. In doing so it must be kept in mind that both these authors have 
neglected not only the contribution of k^^^^ but also the anelasticity and the frequency 

dependence of the Love numbers /c^^^^ Moreover, Dow includes in his analysis also the 
indirect influences of the oblateness of the Earth for 0\, Qi, M2 and A^2- Obviously, 
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in their analyses the LAGEOS II is not present since it was launched only in 1992. 
Another important factor to be considered is the actual sensitivity in measurements of fl 
and uj, in the sense that the eventual discrepancies between the present results and the 
other ones must be not smaller than the experimental error in the Keplerian elements 
if one wish to check the theoretical assumptions behind the different models adopted. 
Carpino has analyzed the inclination and the node of LAGEOS only. His value for the 
important zonal 18.6-year tide is —1087.24 mas, while Table I gives —1,079.38 mas; 
the difference amounts to 7.86 mas, the 0.72 % of the "elastic", frequency-independent 
Carpino's value. Considering that for the node fl of LAGEOS the present accuracy is of 
the order of the mas, 7.86 mas could be in principle detected, allowing for a discrimination 
between the different models adopted in the calculation. For the Ki tide, one of the most 
powerful constituent in perturbing the satellites' orbits, Table I quotes 1,744.38 mas 
against 2144.46 mas of Carpino's result; the gap is 400.08 mas, the 18.6 % of Carpino. 
In the sectorial band, the present analysis quotes for the K2 —92.37 mas and Carpino 
—97.54 mas; there is a difference of 5.17 mas, the 5.3 % of the Carpino's value. It must 
pointed out that there are other tidal lines for which the difference falls below the mas 
level, as is the case for the 9.3-ycar tide. As it could be expected, the major differences 
between the present "anelastic" , frequency-dependent calculations and tli(^ other ones 
based on a single, real value for the Love number ^2 lie in the diurnal band: in it the 
contribution of anelasticity is not particularly relevant, but, as already pointed out, the 
elastic part of k^2i is strongly dependent on frequencies of the tidal spectrum. May be 
interesting to notice that when the calculation have been repeated with the same value 
k2 = 0.317 adopted by Carpino, his results have been obtained again. 

Another feature which characterizes this study is the calculation of the phase lag of 
the solid Earth' s response with respect to the tide generating potential. In Table I, Table 
II and Table III ta.n Simf is always negative and does not exceed 10~^ in absolute value 
showing that the response of the Earth is slightly retarded with respect to the forcing 
lunisolar tidal potential. In absolute value, tan(52o/ is greater than tan 621 f and tan(522/ 
of one order of magnitude, confirming that the behavior of the solid Earth's response is 
more influenced by the anelasticity in the zonal band that in the diurnal and semidiurnal 
ones. 

Up to now the effects induced by the Earth's flattening and the Earth's rotation have 
been neglected. If they have to be analyzed it is necessary to take in eq. (63) the part 



and work out it in the same manner as done for the spherical non-rotating Earth's 

contribution. In applying the transformation to the orbital elements given by cq. (73) 
one has to substitute everywhere I with 1 + 2. So the equations for the perturbations on 
the node and the perigee become 



(93) 
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(94) Xk^mYY. £ ^-^'+2 pqSm-ffl+2 mpq, 

p=0 q=-oo 



p 



(95) X 



^ ^ dGi+2 pq cosidFi+2 

p=0 g 



tl+2mp^^ ^ ^'+' 



sin7/;+2 mpg- 



The corrections induced by the Earth's flattening and rotation, due to the smallness of 
fcj^, have been calculated only for those tidal lines which turned out to be the most 
effective in perturbing the node and the perigee of LAGEOS and LAGEOS II, i.e. the 
zonal 18.6-year tide and the Kx. Adopting the values quoted in ref.[61] for k'l^ and eqs. 
(94)-(95) it is possible to obtain the results summarized in Table IV. The quoted values. 

Table IV. - Corrections, in mas, to the Earth solid tidal perturbations on the nodes of LAGEOS 
and the perigee of LAGEOS II due to the Earth's flattening and Earth's rotation for 1 = 2, p = 
2, g = 0. The values adopted for fc^ are those quoted in ref.[61J. 



Tide 


AfiLAGEOS 


AfiLAGBOS II 


AWLAGEOS II 




055.565 


0.43 


1.28 


0.19 


-0.00094 


165.555 Ki 


-0.97 


-9.49 


0.85 


-0.00074 



with the exception of the node of LAGEOS II, fall below the mas level turning out, at 
the present, undetectable. 

3'5. The mismodelling in the nodes and the perigee. - Concerning the errors in the 
values released in Table I, Table II, and Table III, the major source of uncertainty in 
the perturbative amplitudes lies in the Love numbers A;2^(/). and the orbital injection 
errors Si affecting the inclination i. 



About the latter, by assuming 6i — 0.5 mas [68], wc have calculated 



dA{n) 
di 



Si and 



dAjui) 
di 



Si for 055.565, Ki, and 52 which arc the most powerful tidal constituents in 

perturbing LAGEOS and LAGEOS II orbits. The results are of the order of 10~^ mas, 
so that we can neglect the effect of uncertainties in the inclination determination. 

Concerning the Love numbers fc2rn(/)i from a preliminary point of view, theoreticians 
claim that they are able to compute them with an accuracy of lO"** using the specified 
Earth model. But the Earth models may differ. Probably uncertainty in Earth models 
may bring errors in computation of the Love numbers of the second order about 10~^ ( L. 
Petrov, private communication, 1999). First of all, this means that the level of accuracy 
of 1 mas, i.e. 10~^ rad, on Q and w given by eq. (85) and eq. (86) can be well reached, 
since the uncertainty in Hp{f) amounts to 10^^ m. Second, assuming a mean value of 
0.3 for ^2, it can be stated that the uncertainty in the solid Earth tidal perturbation due 
to a single tide generating potential's constituent amount to almost 0.3 % for O and w. 
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More specifically, we have assessed the uncertainties in them by calculating for cer- 
tain tidal constituents the factor ^^2/^2; ^2 is the average on the values released by the 
most reliable models and 5k2 is its standard deviation. According to the recommen- 
dations of the Working Group of Theoretical Tidal Model of the Earth Tide Commis- 
sion (http://www.astro.oma.be/Dl/EARTH_TIDES/wgtide.html), in the diur- 
nal band wc have chosen the values released in refs. [57, 56] and the two sets of rcf. [61] . For 
the zonal and sectorial bands we have included also the results of ref. [60] . The uncertain- 
ties calculated in the Love numbers span from 0.5% to 1.5% for the tides of interest. 
However, it miist be noticed that the worst known Love numbers are those related to 
the zonal band of the tidal spectrum due to the uncertainties in the anelasticity of the 
Earth' s mantle. These results have been obtained in order to calculate the mismodelled 
amplitudes of the solid tidal perturbations Sfl^^ and Suj^^ for the nodes of LAGEOS 
and LAGEOS II and the perigee of LAGEOS II; they have been subsequently compared 
with the gravitomagnetic precessions over 4 years AfiL^p = 124 mas, AO^t = 126 mas 
and Awlt = —228 mas. In Table V we have quoted only those tidal lines whose mismod- 
elled perturbative amplitudes are greater than 1% of the gravitomagnetic perturbations. 
It turns out that only 055.565 18.6-year and Ki exceed this cutoff. 

Table V. - Mismodeled solid tidal perturbations on the nodes Q of LAGEOS and LAGEOS II 
and the perigee lj of LAGEOS II compared to their gravitomagnetic precessions over 4 years 
Af^Lx = 124 mas Af^Lrj, = 126 mas Awlt = —228 mas. Only those constituents which exceed 
the 1 % cutoff have been quoted. The ratios are in percent and the 5X in mas. 



Tide 


Sk2 
k2 








an" 






055.565 


1.5 


-16.5 


13.3 


30.3 


24 


-21 


9.2 


165.555 Ki 


0.5 


9 


7.2 


-2 


1.6 


10.2 


4.4 



3'6. Gravitational potential of the ocean Earth tides. - One of the early attempts to 
explain the phenomenon of fluid tides [69] in its globality is the equilibrium theory [70, 20] 
by Newton. It is based on the assumption that at every time the free water surface of 
the oceans coincides exactly with the spheroidal equipotential surface at r = ii due 
to the combined action of the Earth's proper gravity and the hmisolar tidal potential. 
Differently stated, the oceanic tidal bulge at every instant coincides exactly with the 
envelope of the forces produced by the tide generating potential at the Earth's surface. 
At every tidal constituents $(/) corresponds an equilibrium partial tide ry/. Referring 
to eq. (61) it can be write 

(96) r]f = = Ai^HpPl^ism cj>) cos{at + mX + ci^). 

For 1 = 2 and m = eq. (96) represents a standing wave which crests at < = 0, while for 
I = 2 and m = 1, 2 it refers to running waves around the oceans from East to West. If 
the equilibrium theory had been valid, it would mean that the water masses acted upon 
by the tidal potential are deprived almost entirely of their inertia in the sense that they 
would adapt instantaneously to their equilibrium positions on the equipotential surface 
which, in fact, changes in time at the frequencies of the lunisolar tidal potential. The 
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reality is quite different, because of tlie complex hydro dynamical behavior of the oceans; 
at every equilibrium partial tide rif corresponds an effective partial tide [71] 

(97) Cf = ^f{d,X)cos[at-Sf{^,X)], 

where £_f{{),\) and (5/(19, A) [72] arc defined as harmonic constants: ^/(i9, A) is half the 
difference between high water and low water for a given place while 6f{'d, A) is the phase 
lag of the ocean tide with respect to the equilibrium one when the phase of the latter 
is calculated at Greenwich meridian A = 0: the delay at any other meridian is obtained 
simply adding mX to df{'&,X). Moreover, Q at every time t crests on the "cotidal" line 
at = 6f{'&,X); in general, referring to a given point (i?, A), e.g., on the shoreline, the 
cotidal lines related to two conscciitivc instants may depart from the shore into the sea 
or vice versa, giving rise to outcoming and incoming waves with respect to that place. 

Such departures from the equilibrium theory are due to several reasons. First of all, 
the fluid elements can freely flow and once they had been put in motion, their notable 
inertia prevent them to change instantaneously their state of motion making them stop 
at the equipotential surface. Second, dissipative phenomena and non-linear interactions 
among the various partial tides and other ocean currents must be taken in account. This 
implies that the response of the water masses to the forcing lunisolar potential presents a 
phase lag with respect to it: the sea surface is not an equipotential one so that tangential 
forces acts on the water particles generating tidal currents. The equilibrium theory may 
be considered an almost good approximation to the reality only at very low frequencies 
of the tidal potential, like as the Moon and the Sun would remain fix in the space with 
respect to the Earth. 

The equations of motion governing the complex hydrodynamics of the tidal currents 
are the so called Lagrange Tidal Equations (LTE) [73, 74, 75, 76, 77], which, in general, 
are not linear. The simplest form for them is obtained considering a spherical, rotating 
Earth entirely covered by an ocean made of a perfect and incompressible fluid; if only 
the tangential motion is considered, neglecting the squares and the cross products of the 
velocity field's components, it is possible to obtain the following linearized LTE 

dC _ 1 Mhvx) d{hsmi& v^) 
^ ^ dt Rsin^^ dX ^ J' 

in which 

1? terrestrial colatitude. 

h depth of the ocean, h = h{'d, A), m. 
C non-equilibrium, ocean partial tide, m. 

and v\ tangential components of the water's velocity field. 
The LTE are linear, and so each ocean partial tide Q can be treated independently of the 
other ones. The fundamental problem of the ocean tides consists of the determination of 
the harmonic constants for every place on the Earth: this can be obtained, in principle, 
resolving the LTE for a given constituent, but this task is in many cases prohibitive due 
to the great calculational complexity needed to obtain realistic results. Another strategy 
consists of employing geodetic artificial satellites, as is done in the realization of the 
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various Earth Gravity Model of the Goddard Space Fhght Center among which EGM96 
[19] is the most recent, or in altimetric measurements [78]. 

In order to calculate the gravitational effect of the water masses raised by the tidal 
forces of a given constituent it can be considered a spherical layer with radius R 

endowed with a bidimensional surface mass density A) = pQ{'&,X), where p is the 
volumetric ocean water density assumed constant, and calculate its potential as 



(101) 



Uf = G 



M/(^,A) 



dT, =G 



I 



In eq. (101) the geodetic convention for the potential is used [1]: in it U is the potential 
usually defined in physics with the sign reversed. Let us rewrite C/ as 

(102) (f = cos{at — 6f) = cos (5/ cos at + ^/ sin 6f sin at. 

The following step consists of expanding eq. (102) in spherical harmonics [35, 36] 

oo / 

(103) 
(104) 



^/ cos6f = X! X! ('^irnf cosmA + bimf sinmA)P™(sin0), 

;=0 rn=0 

sin Sf = X! X! cos'^A + dlmf sinroA)P™(sin^). 



/=0 m=0 



The coefhcients of the harmonic expansion, which have the dimension of lengths, can be 
calculated if the harmonic constants are known by means of the following formulas 

X J ^/cos^/cosTOAP;'"(sin^)dS, 

X y ^/ cos Sf sin mAP/"(sin <j))dJ^, 

X J sinSf cosmXP™{sm(l))dT,, 

^/ sin Sf sinmAP™(sin (j))dYi. 



(105) 


aimf 


(106) 


blmf 


(107) 


Clmf 


(108) 


dlmf 



{2l + l){2~ 


dom) {1 - m)l 


An 


{l + m)\ 


(2Z + 1)(2- 


Som) (l -m)l 


4tt 


{I + m)l 


{2l + l){2- 


Som) {I - my. 


An 


(Z + m)! 


(2Z + 1)(2- 


5om) {l-m)\ 




{l + m)\ 



Eqs. (105)- (108) arc the connection between the harmonic constants obtained through 
oceanographic models and the data acquired by geodetic satellites, in the sense that if 
the coefficients aimf, hmf, cimf, dimf, for a given tidal line, are known from orbital 
analysis, by means of eq. (105)-cq. (108), the essential features of that partial tide may 
be recovered. With eqs. (103)-(104) eq. (102) becomes 



(109) 

in which 

(110) 
(111) 



C/ = X! X! (Cim/ cos mX + Slmf sin mX)P"' (sin ( 



;=o m=0 



Clmf = aimf cos at + Clmf sin at, 
Slmf = blmf cos at + dlmf sin at. 
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By inserting cq. (109) in cq. (101) and expanding in spherical harmonics also the term 



1/ 



, it can be obtained 



(112) Uf = 4ttGRpY^ E WTl { ~ ) ^^""^ "^^'""^ + ^""^ sinmA)Pr(sin<^). 

1=0 m=0 V / 

Ineq. (112), using eq. (llO)-eq. (Ill), one can further expand ((7;^/ cos mA+^jm/ sin mA) 
by introducing the prograde (westwards) and retrograde waves [72] 

(113) Cimf cosmA + Simf sinmA = ^hnf ('^'^ ^ "'0 ^ ^imf ("^'^ ^ ^*)- 



In eq. (113) and in the following, expressions like cos (a ± b) must be interpreted 

as cos (a + &) + A~ cos (a — 6). The sign + refers to the progressive (westwards) waves 
while the sign — is for the waves moving from West to East. From eqs. (llO)-(lll) and 
eq. (113) one has 



(114) 


aflm = + Kflm^ 


(115) 


bflm = S'j^lrn ^ ^flm' 


(116) 


Cflm = S^lm + ^flm ' 


(117) 





Prom these equations, by defining 

(118) K^mf=CLf^ose%^, 

(119) Stf = CtfS^r.e%^, 

it is possible to express the coefficients of the expansion of the harmonic constants 
aimf, bimf, cimf, dim} in terms of C^j, ef^j which can be recovered by geodetic 

O'lmf — dlmf 



satellites 






(120) 






(121) 


^Imf 


COSE;";,/ = 


(122) 






(123) 




sin = 



2 

Cimf + hmf 
2 

-hmf + Cimf 



Eq. (106) and eq. (108) show that, for m = 0, it turns out biof = diof = identically; 
this implies through eqs. (120)-(123) that 

(124) C^f = Cjo/' 
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(125) e;"{jj — Ejo^. 

Eqs. (124)- (125) will be used in the next section. It is interesting to notice that eqs. 
(120)-(123), or also eqs. (114)-(117), show that if it were cimf = dimf = 0, the distinction 
between the prograde and retrograde waves would disappear and only one kind of waves 
would remain, as in the equilibrium theory. Indeed, eqs. (107)-(108) point out that 
Qm/ = dimf = if and only if (5/ = 0, where Sj is the phase lag of the ocean response 
with respect to the lunisolar tidal constituent due to the ocean hydrodynamics. 
Eq. (112), with eq. (113) and eqs. (118)-(119), becomes 

oc I - , „x ;+i (j± 

(126) t//=47rGi?p^^^ - ^cos(at±mA-£,t^)Pr(sin<^). 

Eq. (126) expresses the gravitational potential of the Earth's ocean waters, thought as 
a spherical layer of mass p^/ raised by the action of a given lunisolar tidal constituent of 

circiilar frequency cr. 

But these enormous water masses act upon the sea floor and the whole of the solid 
Earth which are deformed and attracted by them. This effect has also to be considered. 

It can be done according to cq. (60) and introducing the so called load Love numbers fc^ 
[79, 35, 80]. The global response of the Earth oceans to the forcing tidal constituent of 
frequency a can be written as 

X 



=0 m=0 + 



r J 21 + 1 



(127) xCtf cos {at ±mX- etf)Pr {sin ct>), 

with Uf < Uoc, i-e. A: < 0. This fact can be intuitively guessed by thinking about U 
in terms of shift of the equipotential surfaces h = U/g: if the sea floor is raised by the 
tidally driven water masses the total shift is smaller than it would be if we considered 
the solid Earth as completely rigid. 

3T. Ocean tidal orbital perturbations. - Eq. (127) can be fruitfully rewritten in a 

form more suitable for calculating the perturbations induced on the Kcplcrian orbital 
elements of artificial satellites [81, 82]. To this aim the following quantities are defined 

(128) Atf=4nGRp^l^Ctf, 

(129) C/=Ai/Cos(T^t±e±„^), 

(130) = sin {Tat ± e^f). 

Uf becomes 

(131) C^/ = E E E (7) ^ (C,^^cosmA + <S,^^sinmA)Pr(sin<^). 
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Eq. (131) is formally equal to the expression of the static gravitational potential of the 
Earth worked out by Kaula. So, with the same mathematical menagerie used for the 
geopotcntial in ref.[l], it is possible to write cq. (131) in terms of the Keplerian orbital 
elements of a test body in the field of the Earth as 



(132) Uf 



OO I — y p\ l + l I +00 

;=0 m=0 + ^ p=0 g=-oo 



in which 
(133) 



^fimpg = (1- 2P)'^ +{l-2p + q)M + m(n - 9) ± {at - ef^^). 



The equations for the tidal ocean perturbations may be worked out as already done for 
the solid Earth tides in subsect. 4. At first order, one obtains 
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It should be noticed that the frequencies of the perturbations given by eq. (136) are 
different, in general, from the frequencies of the solid Earth tidal perturbations given 
by eq. (88). While for the solid tides the diurnal modulation due to 9 cancels out 
automatically if one considers those terms in which ji = m, for the ocean tides, in 
general, this docs not happen because of the presence of the Eastwards waves due to the 
non-equilibrium pattern of the ocean tidal bulge. 

Long periodicities can be obtained considering those combinations of I, p and q for 
which I — 2p + q = holds and retaining the Westward progradc terms with ji = m. 
Only in this way in eq. (136) the contributions of 9 are equal and opposite, and cancel 
out. With these bounds on I, m, p and q the frequencies of the perturbations of interest 
become 



(137) 



/p = f / + (/ - 2p)oj + mQ. 
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It is worthwhile noting that the frequencies fp are identical to those of solid tidal per- 
turbations [35, 51], so that the satellites cannot distinguish one effect from the other, a 
feature which will turn out to be important in the recovery of the coefficients C;^„y. 

For Z = 2, as for the solid tides, l — 2p = holds. The frequencies of the perturbations 
are, in this case, given by 

(138) fp = tf + mh. 

If I is odd the situation is different because now all the terms with p running from 
to I must be considered and q is different from zero: q = 2p — I. So, the rates of the 
perturbations include the contribute of io. 

Up to now the prograde wave terms only have been considered in order to deal with 
long period perturbations. Yet, there is a case in which also the retrograde terms fit into 
this scheme; it is possible to show that for I even and m = the frequencies of prograde 
and retrograde terms coincide, in absolute values, and the perturbation amplitudes are 
twice the amplitudes of the prograde terms only. Indeed, according to eqs. (134)-(135), 
Af2/, Aw/ are proportional to 



(139) E4™T 



+ 

in this case / — m is even and retaining only those terms for which I — 2p = one has 

(140) 7/1oio = ±('^*-4o)' 

(141) fp = ±a. 
Expanding eq. (139) it can be obtained 

^/(O^ sin {at - 4;o) + ^JioZ^ sin [-{at - eJ^q)] = 



(142) = A+Q^ sin {at - e+J + ^/(o^ sin i^rt -e^J. 

Since A'^^Q = Aj^g and Ej^q = eJ^q, as shown by eq. (124) and eq. (125), eq. (142) 
becomes 



(143) _i^£sin(at-4o). 



Eqs. (134)-(135) have been adopted in order to compute the amplitudes of the ocean 

tidal perturbations on the Keplerian elements Vl for LAGEOS and LAGEOS II and u> 
for LAGEOS II which pertain directly the Lense-Thirring experiment. Also in this 
section the perturbations of second order [65] due to the oblateness of the Earth have 
not been considered, as already done for the solid tides. The inclination and eccentricity 
functions used are those quoted in [1] . For the numerical values of the various geophysical 
parameters which figure in eqs. (134)-(135) the lERS standards [56] have been used, while 
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the EGM96 gravity model [19] has been adopted for the choice of the tidal constituents 
and their coefBcients in the ocean tidal potential Uf. The calculation have been 
performed, considering only the progressive waves, for the following tidal lines 

• Mm (065.455), Sa (056.554), M/ (075.555), Ssa (057.555); 

• Ki (165.555), Oi (145.555), Pi (163.555), Qx (135.655); 

• if 2 (275.555), Af2 (255.555), S2 (273.555), N2 (245.655), T2 (272.556). 
For each of these tidal constituents the following terms have been calculated 

1 = 2, p = 1, q = because the eccentricity functions Gipq{e) for = 0, q = —2 and 
p = 2, q ~ 2 vanish. 

1 = 2), p = I, q = —1 and / = 3, p = 2, q = I because G30-3 and G333 are not quoted 
in [1] due to their smallness: indeed, the Gipq{e) are proportional to el^L 
1 = 4:, p = 2, q = because the other admissible combinations of I, p and q give rise to 
negligible eccentricity functions. So, also for ^ = 4 the condition ^ — 2p = 0, in practice, 
holds and the constituents of degree I = 2 and I = 4 generate detectable perturbations 
with identical periods. 

Similar analysis can be found in [67, 35] . It should be noticed that when Christodoulidis 
performed his study, which is relative to only five constituents, neither the LAGEOS nor 
the LAGEOS II were in orbit, while Dow has sampled the tidal spectrum for LAGEOS in 
a poorer manner with respect to this study in the sense that, for each constituent, only 
the terms of degree / = 2 have been considered with the exception of the Ki whose I = 4 
contribution has been also analyzed. Moreover, when these works have been realized 
there were a few coefficients Cj^^ available with relevant associated errors. 

Some explanations are needed about the determination of the coefficients C^j:. In 
EGM96 the geopotential is recovered through both altimeter and surface gravity infor- 
mation, and satellites data. The pool of near-Earth satellites employed, in general, are 
not perceptibly perturbed by the entire tidal spectrum, but they are sensitive only to 
some certain tidal lines, depending on the features of their orbits. Moreover, on a large 
enough temporal scale they cannot distinguish between solid and prograde ocean tidal 
perturbations because their frequencies are the same. Finally, there are also other pe- 
riodic physical phenomena different from the tides that affect the orbits of the gcxxletic 
satellites; in many cases their periodicities are similar to that of the tides, particularly 
in the zonal band. These facts imply that in the analytical expressions of the perturba- 
tions it is necessary to assume as known a priori from various reliable models some solid 
or ocean tidal terms and consider variable the other ones in order to adjust them by 
means of the experimentally determined values of the perturbations; all the constituents 
considered, both those held fixed and the other recovered, must be capable to influence 
perceptibly the satellites employed. The strategy followed in EGM96 has consisted in 
adopting the frequency-dependent Wahr model for the solid tides [52] with its values for 
the Love numbers and the Hp{f) of Cartwright and Edden, calculating a certain number 
of ocean tidal terms by means of oceanographic models, when it has been possible to solve 
the LTE, or by means of other algorithms [37], considering also the retrograde waves, 
and choosing the 13 ocean tidal terms listed above to be adjusted. The values recovered 
for the coefficients C^Qf account for also the retrograde terms because, as previously 
noticed, their periods are equal to those of prograde ones. The terms whose value has 
been considered given a priori constitiite the so called background; it necessarily contains 
terms up to sixth degree because the attenuation due to altitude makes the near-Earth 
satellites almost insensitive to larger degree terms: indeed, in eqs. (134)-(135) the 
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are proportional to 
(144) 



rV+' 1 



21 + 1 



The values obtained for the coefficients C'j^j: are, in general, biased by the effects of the 
anelasticity of the solid Earth's mantle and by all other phenomena which have not been 
explicitly modeled in the background. For example, in the values rc;covered for the Sa are 
included climatological effects which have not gravitational origin; in the 5*2 coefficients 
are also included the variations of the atmospheric pressure due to the atmospheric tides. 

3'8. Discussion of the numerical results. In Table VI, Table VII and Table VIII 
the present results for the nodes Q of LAGEOS and LAGEOS II, and the argument of 
perigee w for LAGEOS II are quoted. 

Table VI. - Perturbative amplitudes on the node f2 of LAGEOS due to ocean tides. P indicates 
the periods in days, A the amplitudes in mas and E the percent error in the C^^j: . The values 
employed for them and the related errors are those quoted in EGM96 model. 



Tide 


1 = 2, p = 


1, 9 = 




I = 3, p 


= !,<! = 


-1 


I = 3, p 


= 2, 9 = 1 




! = 4, p = 


2, 9 = 


P 


A 


B 


P 


A 


E 


P 


A 


B 


P 
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Prom an accurate inspection of Table VI, Table VII it is possible to notice that, for the 

nodes 17, only the even degree terms give an appreciable contribute; the I = 3 terms are 
totally negligible. This so because All/ is proportional to the Gipq{e) functions which, 
in turn, are, in general, proportional to el^l; in this case the eccentricity functions are 
G(e)3i_i = e(l - e^)^^/^ = G(e)32i and the eccentricities of LAGEOSsatellites are of 
the order of 10~-^ or less. 

Among the long period zonal tides, the Solar annual tide Sa (056.554) exerts the most 
relevant action on the nodes, with an associated percent error in the amplitudes of 6.7 
%. It is interesting to compare for this tidal line the ocean tidal perturbations of degree 
I = 2 A£^gEOs(^) = -20.55 mas, A^^^^qs ii(^) = 37.7 mas, with those due to the 
solid Earth tides Alageos(^) = 9.96 mas, A^°ageos ii(^) = -18-28 mas. The ocean 
amplitudes amount to 204 % of the solid ones, while for the other zonal constituents they 
vary from 9.9 % for the Mm to the 18.5 % of the Ssa- This seems to point toward that 
the recovered value of for the Sa is biased by other climatological effects than the 
tides; indeed, for all the other tidal lines, zonal or not, the ocean tidal perturbations of 
degree / = 2 amount to almost 10 % of the solid Earth tides perturbations. 
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Table VII. - Perturbative amplitudes on the node of LAGEOS II due to ocean tides. P 
indicates the periods in days, A the amplitudes in mas and E the percent error in the Cp^^ ■ The 
values employed for them and the related errors are those quoted in EGM96 model. 
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Table VIII. - Perturbative amplitudes on the perigee u> of LAGEOS II due to ocean tides. P 
indicates the periods in days, A the amplitudes in mas and E the percent error in the Cp^^. The 
values employed for them and the related errors are those quoted in EGM96 model. 
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The terms of degree I = 2 of the tesseral tides Ki (165.555) and Oi (145.555) induce 
very large perturbations on the node of LAGEOS and, to a lesser extent, of LAGEOS 
II; Table VI and Tabic Vll quote 156.55 mas, 151.02 mas for the former and —35.69 
mas, —34.43 mas for the latter. The associated percent errors are 3.8 % and 2.9 %, 
respectively. For the terms of degree I = 4, which have the same periods of those of 
degree I = 2, the situation is reversed: Table VI and Table VII quote 4.63 mas and 3.54 
mas for LAGEOS and 41.58 and 31.81 mas for LAGEOS II. The percent errors associated 
with Ki and Oi, for I = 4, are 3.9 % and 5.7 % respectively. For all the tesseral lines 
investigated the ocean tidal perturbations of degree I = 2 are, in general, the 10 % of the 
solid tidal perturbations. 

Among the sectorial tides, the most relevant in perturbing the nodes of LAGEOS 
satellites, on large temporal scales, is the M2 (255.555): Table VI and Table VII quote 
-75.59 mas {I = 2) and -77.40 mas {I = 4) for LAGEOS, and -75.65 mas {I = 2) and 
-48.07 mas {I = 4) for LAGEOS 11. The associated percent errors are 0.9 % for Z = 2 
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and 42.8 % for i = 4. The M2 ocean perturbation of degree / = 2 amounts to 10.3 % of 
the corresponding sohd tide. The same holds for the other sectorial tides of degree I = 2. 

About the periodicities of these perturbations in relation to the detection of Lense- 
Thirring effect, the same considerations already exposed for the solid tides hold also 
in this case. Moreover, it must be pointed out the aliasing role played by the zonal 
18.6-ycar and 9.3-ycar tides in the extracting any secular effect, like the gravitomagnctic 
precession, from a record whose duration is shorter than their very long periods. Recently 
both Starlette and LAGEOS SLR satellites passed their 19th year in orbit and this span 
of time is now adequate to get reliable information about these tides [83, 84] due to 
their slow frequencies they can be correctly modelled in terms of the equilibrium theory 
through the HJ" coefficients and a complex Love number accounting for the anelasticity of 
the mantle. So, concerning them the results quoted for the solid tides can be considered 
adequately representative. In Table VIII the amplitudes of the perturbations on the 
argument of perigee w for LAGEOS II are quoted. For this orbital element the factor 



to which Aw j is proportional makes the contributions of the I = 3 terms not negligible. 
For the even degree terms the situation is quite similar to that of f2 in the sense that the 

most influent tidal lines arc the Sa^ Ki, Oi and Af2- 

Once again, among the long period tides the Sa exhibits a characteristic behavior. Indeed, 
its Z = 3 contributions are much stronger than those of the other zonal tides. This fact 
could be connected to the large values obtained in its C^qj coefficient and people believe 
that it partially represents north to south hemisphere mass transport effects with an 
annual cycle nontidal in origin. The I = 3 terms present, in general, for the perigee of 
LAGEOS II a very interesting spectrum also for the tesseral and sectorial bands: for 
LAGEOS II there are lots of tidal lines which induce, on large temporal scales, very 
relevant perturbations on u, with periods of the order of an year or more. In particular, 
it must be quoted the effect of Ki line for p = 1, q = —I: the perturbation induced 
amount to -1,136 mas with period of -1,851.9 days. These values arc comparable to the 
effects induced by the solid Earth tidal constituent of degree I — 2 on the node of 
LAGEOS. Comparing the degree I = 2 terms with those of the solid tides, it can be 
noticed that also for the perigee the proportions are the same already seen for the nodes. 

3'9. The mismodelling in the nodes and the perigee. - Concerning the mismodelling 
on the ocean tidal perturbations, the main source of uncertainties in them is represented 
by the load Love numbers fcj and the coefficients C^^. In regard to the ocean loading, 

the first calculations of the load Love numbers can be found in refs.[79]. Pagiatakis, in 
ref.[80], in a first step has recalculated fc; for an elastic, isotropic and non-rotating Earth: 
for I < 800 he claims that his estimates differ from those by Farrell in ref.[79], calculated 
with the same hypotheses, of less than 1%. Subsequently, he added to the equations, 
one at a time, the effects of anisotropy, rotation and dissipation; for low values of I their 
eff'ects on the results of the calculations amount to less than 1%. It has been decided 
to calculate ^^p^ 6k', of the perigee of LAGEOS 11 for Ki I = 3 p = 1 which turns 

out to be the most powerful ocean tidal constituent acting upon this orbital element. 
First, we have calculated mean and standard deviation of the values for fcg by Farrell 
and Pagiatakis obtaining 5k^/k^ = 0.9%, in according to the estimates by Pagiatakis. 
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Then, by assuming in a pessimistic way that the global effect of the departures from 
these symmetric models yield to a total Sk,^/k,^ = 2%, we have obtained 5ui^^ = 5.5 
mas which corresponds to 2% of AciJlt over 4 years. Subsequently, for this constituent 
and for all other tidal lines we have calculated the effect of the mismodelling of 
as quoted in EGM96 [19]. In Table IX we compare the so obtained mismodelled ocean 
tidal perturbations to those generated over 4 years by the Lense-Thirring effect. It turns 

Table IX. - Mismodeled ocean tidal perturbations on the nodes O, of LAGEOS and LAGEOS 
II and the perigee cu of LAGEOS II compared to their gravitornagnetic precessions over 4 years 
AOlt = 124 mas A^lt = 126 mas Awlx = —228 mas. The effect of the ocean loading has been 
neglected. When the 1% cutoff has not been reached a - has been inserted. The values quoted for 
Ki I = 3 p = 1 includes also the mismodelling in the ocean loading coefficient kg assumed equal 
to 2%. The ratios are m percent and the mismodelled amplitudes 5X in mas. 
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18 


7.9 


Ki 1=4 p=2 q=0 


3.9 






1.6 


1.2 






Pi 1^3 p=l q^-1 


18.5 










5.3 


2.3 


Pi 1=3 p=2 q=l 


18.5 










6.4 


2.8 


K2 1=3 p=l q=-l 


5.5 










11.7 


5 


K2 1=3 p=2 q=l 


5.5 










4.8 


2 


S2 1=3 p=l q=-l 


7.1 










6.9 


3 


S2 1=3 p=2 q=l 


7.1 










4.4 


1.9 


T2 1=3 p=l q=-l 


50 










2.6 


1.1 



out that the perigee of LAGEOS II is more sensitive to the mismodelling of the ocean 
part of the Earth response to the tide generating potential. In particular, the effect of 



Ki I = 3 p = 1 q = —1 is relevant with a total 5uj = 




mas amounting to 28.3 % of Ac^lt over 4 years. 

3T0. Conclusions. - Concerning the orbital tidal perturbations on LAGEOS and 
LAGEOS II the following improvements with respect to the previous works have been 

reached 

• The node and the perigee have been considered 

• The analysis has been extended also to LAGEOS II, laimched in 1992 

• Concerning the solid tides, the most recent available frequency-dependent Love 
numbers k2 have been used instead of a single-valued Love number ^2 = 0.317. Moreover, 
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the latitude-dependence of the Love number fc+ has been considered for some selected 
tidal lines 

• About the ocean tides, their orbital perturbations have been extensively calculated 
for / = 2,3,4 by using the most recent available model EGM96. Interesting are the 
results obtained for the Ki, l = 3m = lp=lq = —l oceanic constituent. Indeed it 
induces on the perigee of LAGEOS II perturbations whose amplitudes are of the order of 
thousands of mas and the periods amounts to some years as the solid tidal perturbations 

• An evaluation of the mismodelling on such orbital perturbations, for certain tidal 
constituents, have been performed 

The calculations performed here have a general validity because they can be extended 
to any artificial satellite. The results presented here for the LAGEOS can be used, from 
one hand, in improving the modelling of their orbital perturbations and, from another 
hand, for the correct evaluation of the error budget in any space-based experiment de- 
voted to the measurement of some particular feature of the Earth' s gravitational field. 

In the context of the general relativistic Lense-Thirring experiment, the calculations 
performed here have been used, in view of a refinement of the error budget, in order 
to check preliminarily which tidal constituents are really important in perturbing the 
combined residuals so to fit and remove them from the data, if possible, or, at least, to 
evaluate the systematic error induced by them. Table eq. (V) and Table eq. (IX) show 
that, over a 4 years time span the nodes of the two LAGEOS are sensitive to the even 
components of the 18.6-year line, the Ki and the Sa at a 1% level at least. Moreover, 
the perigee of LAGEOS II turns out to be very sensitive to the / = 3 part of the ocean 
tidal spectrum. 

The results obtained in this section for the nodes of LAGEOS and LAGEOS II and 
the perigee of LAGEOS II will be the starting point for the numerical simulations with 
MATLAB, described in section 4 and section 6, performed to assess quantitatively the 
direct impact of the tidal perturbations on the gravitomagnetic Lense-Thirring experi- 
ment, currently implemented, and on the proposed gravitoelectric experiment. 

4. — The impact of the Earth tides on the determination of the Lense-Thirring 

effect 

4"1. Introduction. - According to [14], the Lense-Thirring effect could be detected, 
in the field of the Earth, by analyzing the orbits of the two laser-ranged LAGEOS and 
LAGEOS II satellites. The observable adopted is the combination of orbital residuals 
of eq. (5): the determination of the scaling parameter /xlt, 1 in General Relativity 
and in Newtonian mechanics, is influenced by a great number of gravitational and 
non-gravitational perturbations acting upon LAGEOS and LAGEOS II. Among the per- 
turbations of gravitational origin a primary role is played by the solid Earth and ocean 
tides. Their effects on the orbital elements of LAGEOS and LAGEOS II has been ex- 
tensively analyzed in section 3. 

In [14] it is claimed the combined residuals y = S^l^^ + ci(5rig^p + C2SijJ^^p allow one 
to cancel out the static and dynamical perturbations of degree I = 2,4: and order m = 
of the terrestrial gravitational field; however, this is not so for the tesseral (to = 1) and 
sectorial (to = 2) tides. This section aims to assess quantitatively how the solid and 
ocean Earth tides of order to = 0, 1, 2 affect the recovery of /xlt in view of a refinement 
of the error budget of the gravitomagnetic LAGEOS experiment. 

Concerning the zonal tides, the results obtained in section 3 for the amplitudes of their 
perturbations have been used directly in eq. (5) in order to test in a preliminary way if 
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the I = 2,4 m = tidal perturbations really cancel out. In regard to the tesseral and 
sectorial tides the analysis of their impact on the Lense-Thirring measurement has been 
done by simulating the real residual curve and analyzing it. The Tobs chosen span ranges 
from 4 years to 8 years. This is so because 4 years is the length of the latest time series 
actually analyzed in [15] and 8 years is the maximum length obtainable today (Summer 
2000) because LAGEOS II has been launched in 1992. The analysis includes also the 
long-period signals due to solar radiation pressure and the J3 geopotential zonal harmonic 
acting on the perigee of LAGEOS II. In our case we have a signal built up with the secular 
Lense-Thirring trend('') and a certain number of long-period harmonics, i.e. the tesseral 
and sectorial tidal perturbations and the other signals with known periodicities. 

The part of interest for us is the secular trend while the harmonic part represents the 
noise. We address the problem of how the harmonics affect the recovery of the secular 
trend on given time spans Tobs and various samplings At. 

Among the long-period perturbations we distinguish between those signals whose pe- 
riods are shorter than Tobs and those signals with periods longer than Tobs- While the 
former average out if Tobs is an integer multiple of their periods, the action of the latter 
is more subtle since they could resemble a trend over temporal intervals too short with 
respect to their periods. They must be considered as biases on the Lense-Thirring deter- 
mination affecting its recovery by means of their mismodelling. Thus, it is of the utmost 
importance that we reliably assess their effect on the determination of the trend of the 
gravitomagnetic effect. It would be useful to direct the efforts of the community (geo- 
physicists, astronomers and space geodesists) towards the improvement of our knowledge 
of those tidal constituents to which /xlt turns out to be particularly sensitive (for the 
LAGEOS orbits). 

This investigation will quantify unambiguously what one means with statements like: 
A/itidos < X%iJ,ijT- In this section we shall try to put forward a simple and meaningful 
approach. It must be pointed out though that it is not a straightforward application of 
any exact or rigorously proven method; on the contrary, it is, at a certain level, heuristic 
and intuitive, but it has the merit of yielding reasonable and simple answers and allowing 
for their critical discussion. 

The section is organized as follows. In subsect. 2 the effects of the / = 2, 4 m = 
tidal constituents on the combined residuals is examined; it turns out that, not only they 
affect it at a level << 1%, but this feature also extends, within certain limits, to the 
/ = 3, m = part of the tidal spectrum. In subsect. 3 the features of the simulation 
procedure of the observable curve are outlined. In subsect. 4 the effects of the harmonics 
of the order m = 1 , 2 with period shorter than 4 years are examined by comparing the 
least squares fitted values of /ilt in two different scenarios: in the first one the simulated 
curve is complete and the fitting model contains all the most relevant signals plus a 
straight line, while in the second one some selected harmonics are removed from the 
simulated curve which is fitted by means of the straight line only. Subsect. 5 addresses 
the topic of the harmonics with periods longer than 4 years. Subsect. 6 is devoted to 
the conclusions. 



C) In fact, the combined residuals are affected also by the secular contribution due to i = 6, 8, .. 
zonal terms of the geopotential. Their effect amounts to almost 13% of the Lense-Thirring trend 
[15]. See also Appendix A. 
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4'2. Systematic error due to the zonal tides. - The combined residuals by Ciufolini 
are useful since they should vanish if calculated for the even zonal contributions C20 and 
C40 of the gcopotential [14]. More precisely, the right side of cq. (5) should become equal 
to zero if the left side were calculated for any of these two even zonal contributions, both 
of static and dynamical origin; the nearer to zero is the right side, the smaller is the 
systematic uncertainty in /ilt due to the contribution considered. 

In order to test this important feature for the case of tides, in a very conservative 
way the results obtained in section 3, Table I, Table II, Table III for the solid Earth 
tides and Table VI, Tabic VII, Table VIII for the ocean tides, have been used in cq. (5) 
by assuming, for the sake of clarity and in order to make easier the comparison with 
ref.[66], an observational period of 1 year and the nominal values of the calculated tidal 
perturbative amplitudes, as if the zonal solid and ocean tides were not at all included 
in the GEODYN II dynamical models so that the residuals should account entirely for 
them. 

The results are released in Table X and Table XI. 

Table X. - Contribution of the even zonal solid tidal constituents to A/ilt by means of the 
formula SQ^ + Jfi" x 0.295 - Sd>" x 0.35 = 60.2 x nur for I = 2, m = 0, p = 1, q = 0. 



Tide 


A{Qi) (mas) 


A{i}u) (mas) 


A{uJu) (mas) 




055.565 


-1,079.38 


1,982.16 


-1,375.58 


-0.219 


055.575 


5.23 


-9.61 


6.66 


1.06 X 10"^ 


056.554 Sa 


9.95 


-18.28 


12.69 


2.02 X 10"^ 


057.555 S,a 


31.21 


-57.31 


39.77 


6.33 X 10-3 


065.455 Af„, 


5.28 


-9.71 


6.74 


1.07 X 10-3 


075.")").") Mj 




-9.08 


G.:i 


1 X 10 



Table XI. - Contribution of the even zonal ocean tidal constituents to A/ilt by means of the 
formula 60} + Jfi" x 0.295 - Sdi" x 0.35 = 60.2 x /xlt /or ? = 2, m = 0, p = 1, q = 0. 



Tide 


A{Qi) (mas) 


A{Qu) (mas) 


A(a»ii) (mas) 


A/iLT 


056.554 Sa 


-20.55 


37.71 


-26.17 


-3.68 X 10-3 


057.555 S^a 


-5.98 


10.98 


-7.62 


-1.28 X 10-3 


065.455 Mm 


-0.54 


1 


-0.69 


-8.78 X 10-^ 


075.555 Mf 


-0.62 


1.13 


-0.78 


-1.73 X 10-* 



The values of A/ilt quoted there for the various zonal tidal lines may be considered 
as the systematic error in /xlt due to the chosen constituents, if considered one at a time 
by neglecting any possible reciprocal correlation among the other tidal lines. Table X 
and Table XI show that the percent error in the general relativistic value of /xlt due to 
the 18.6-year tide, the most insidious one in recovering the LT since it superimposes to 
the gravitomagnetic trend over time spans of a few year, amounts to 21.9 %, while for 
all the other zonal tides it decreases to 0.1 % or less. This means that, even if neglected 
in the satellite orbit determination models, the I = 2 m = tides, with the exception of 
the 18.6-year tide, do not affect the recovery of /xlt by means of the combined residuals. 
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It is interesting to compare the present results to those released in ref.[66] for the 
18.6-year tide. The value -0.219 due to the solid component for /xlt quoted in Table X 
must be compared to -0.361 in rcf. [66], with an improvement of 39.3 %. In the cited work 
there is no reference to any estimate of the mismodelling of the 18.6-year tide, so that 
we have used the nominal tidal perturbative amplitudes released in it: A{Q}) = —997 
mas, A(^l^^) = 1, 805 mas and A{uj^^) = — 1, 265 mas. These figures for the perturbative 
amplitudes due to the solid Earth tide of 18.6-year are notably different from those 
quoted in the present study. In ref.[66] the theoretical framework in which they have 
been calculated (F. Vcspc, private communication, 1999) is based on the assumption of 
a spherical, static, elastic Earth with a single nominal value of k2 = 0.317 used for the 
entire tidal spectrum. The inclusion of the tiny corrections due to the Earth's flattening 
and rotation on the perturbative amplitudes of Q and lu for the 18.6-ycar tide could allow 
one to slightly improve the related uncertainty in /zlt; it would amount to 20.6 %. But 
since the present-day accuracy in laser ranging measurements could hardly allow one to 
detect these small effects, their utilization in eq. (5) is debatable. 

Remember that the result quoted for the 18.6-year tide is obtained in the worst 
possible case, i.e. a time span of only 1 year and the assumption that the residuals have 
been built up by neglecting completely the zonal tides in the dynamical models used. If, 
more realistically, we calculate eq. (5) with the mismodclled amplitudes quoted in the 
flrst row of Table V for the 18.6-year tide we obtain, over 1 year, A/i = —3.51 x 10^'^. This 
strongly highlights that many efforts, either theoretical or experimental, must be done in 
order to modelling as accurately as possible such a constituent, so that it can be included 
in the nominal background of the orbital determination softwares like GEODYN II at a 
satisfactory level of accuracy. The calculations performed here point toward this goal in 
the sense that, if we put our values for the perturbative amplitudes due to the 18.6-year 
tide in the models of, e.g., GEODYN II and subsequently build up the orbital residuals, 
we expect that the contribution of such semisecular constituent to A/xlt amounts to the 
value quoted here. 

Even though a cancellation is not expected as for the first two even zonal constituents, 
by calculating the left hand side of eq. (5) for the other tides yields, at least, an order of 
magnitude of their effect on /xlt- An interesting, unpredicted feature stands out for the 
odd zonal ocean tides. The contribution oil — 3 zonal ocean tidal nominal perturbations 
over 1 year to A/xlt can be found in Table XII and Table XIII. 



Table XII. - Contribution of the odd zonal ocean tidal constituents to A/ult by means of the 
formula Stli + 50" x 0.295 - to" x 0.35 = 60.2 x nur for I = 3, m = 0, p = 1, q = -1. 



Tide 


A{Qi) (mas) 


A{Qii) (mas) 


A{u!u) (mas) 


A/iLT 


056.554 Sa 


-0.063 


0.13 


-114.35 


0.66 


057.555 Ssa 


-9 X 10"^ 


0.028 


-22.95 


0.133 


065.455 Mm 


-4 X lO""* 


1 X 10"^ 


-1.53 


-8.93 X IQ-^' 


075.555 Mf 


-1 X lO-* 


7 X 10"* 


-0.58 


3.41 X 10"^ 



By inspecting them it is clear that the sensitivity of perigee of LAGEOS II to the 
I = 3 part of the ocean tidal spectrum may affect the recovery of the Lense-Thirring 
parameter /xlt by means of the proposed combined residuals, especially as far as Sa and 
Ssa ai^e concerned. This fact agrees with the results of Table IX which tells us that 
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Table XIII. - Contribution of the odd zonal ocean tidal constituents to A/^lt by means of the 
formula Stl^ + 50" x 0.295 - Suj" x 0.35 = 60.2 x /xlt /or Z = 3, m = 0, p = 2, q = l. 



Tide A{ni) (mas) A(f2ii) (mas) A{loii) (mas) A/iLx 

056.554 Sa 0.047 -0.36 297.34 -1.72 

057.555 Ssa 7 x lO'^ -0.044 36.07 -0.209 

065.455 M„ 4 x 10"* -2 x 10"^ 1.642 -9.55 x IQ-^ 

075.555 Mf 1.79 x 10""^ -7 x 10""^ -0.60 -3.52 x 10"^ 



the mismodelled parts of Sa and Ssa are not negligible fractions of Awlt. However, if 

the mismodelled amplitudes are employed in eq. (5) it can be seen that, over 1 year, a 
cancellation of the order of 10~^ {Sa I = 3 p = 2) and 10~^ {Sa I = 3 p = 1; Ssa I = 2>p = 
1, 2) takes place. The contributions of the mismodelling on Mm and M/ are completely 
negligible. So, we can conclude that also the I = 3 part of the zonal ocean tidal spectrum 
may not affect the combined residuals in a notable manner if the I — 3 part of Sa and 
Ssa is properly accounted for. 

4'3. The simulated residual signal. - The first step of our strategy was to generate 
with MATLAB a time series which simulates, at the same length and resolution, the real 
residual curve obtainable from y = 6^\^^ + ciJfig^p + c^Suj^J^^ through, e.g., GEODYN. 
This simulated curve ("Input Model" - IM in the following) was constructed with 

• The secular Lcnsc-Thirring trend as predicted by the General Relativity(^) 

• A certain number of sinusoids of the form SAi- cos {"^t + (j)k) with known periods 
Pk, fc = 1, ..N simulating the mismodelled tides and other long-period signals which, to 
the level of assumed mismodelling, affect the combined residuals 

• A noise of given amplitude which simulates the experimental errors in the laser- 
ranged measurements and, depending on the characteristics chosen for it, some other 
physical forces. 

In a nutshell 



(14^ = LT -|- [mismodelled tides] -|- [other mismodelled long period signals] -|- [noise]. 

The harmonics included in the IM are the following 

• Ki, / = 2 solid and ocean; node of LAGEOS (P=l,043.67 days) 

• Ki, 1 = 2 solid and ocean; node and perigee of LAGEOS II (P=-569.21 days) 

• Ki, 1 = 3, p=l ocean; perigee of LAGEOS II (P=-l,851.9 days) 

• Ki, 1 = 3, p^2 ocean; perigee of LAGEOS II (P=-336.28 days) 

• K2, 1 = 3, p=l ocean; perigee of LAGEOS II (P=-435.3 days) 

• K2, 1 = 3, p = 2 ocean; perigee of LAGEOS II (P=-211.4 days) 

• 165.565, I = 2 solid; node of LAGEOS (P=904.77 days) 

• 165.565, I = 2 solid; node and perigee of LAGEOS II (P=-621.22 days) 

• 52, 1 = 2 solid and ocean; node of LAGEOS (P=-280.93 days) 

• ^2, I = 2 solid and ocean; node and perigee of LAGEOS II (P=-111.24 days) 



{^) Remember that in the dynamical models of GEODYN II it was set purposely equal to so 
that the residuals absorbed (contained) entirely the relativistic effect [66] . 
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• S2, 1 = 3, p=l ocean; perigee of LAGEOS II (P=-128.6 days) 

• S2, I = 3, p = 2 ocean; perigee of LAGEOS II (P=-97.9 days) 

• Pi, 1 = 2 solid and ocean; node of LAGEOS (P=-22L35 days) 

• Pi, 1 = 2 solid and ocean; node and perigee of LAGEOS II (P=-138.26 days) 

• Pi, 1 = 3, p=l ocean; perigee of LAGEOS II (P=-166.2 days) 

• Pi, ; = 3, p = 2 ocean; perigee of LAGEOS II (P=-118.35 days) 

• Solar Radiation Pressure, perigee of LAGEOS II (P=-4,241 days) 

• Solar Radiation Pressure, perigee of LAGEOS II (P=657 days) 

• perigee odd zonal C30, perigee of LAGEOS II (P=821.79) 

In the following the signals due to solar radiation pressure will be denoted as SRP(P) 
where P is their period; the effects of the eclipses and Earth penumbra have not been 
accounted for. Many of the periodicities listed above have been actually found in the 
Fourier spectrum of the real residual curve [66]. Concerning Ki I = 3 p = 1 and 
SRP(4,241), see subsect. 5. 

When the real data arc collected they refer to a unique, unrepeatable situation char- 
acterized by certain starting and ending dates for Tobs- This means that each analysis 
which could be carried out in the real world necessarily refers to a given set of initial 
phases and noise for the residual curve corresponding to the chosen observational period; 
if the data are collected over the same Tobs shifted backward or forward in time, in gen- 
eral, such features of the curve will change. We neither know a priori when the next real 
experiment will be performed, nor which will be the set of initial phases and the level of 
experimental errors accounted for by the noise. Moreover, maybe the dynamical models 
of the orbit determination software employed are out of date in regard to the perturba- 
tions acting upon satellites' orbits or do not include some of them at all. Consequently, 
it would be incorrect to work with a single simulated curve, fixed by an arbitrary set of 
SAk, 4>k and noise, because it could refer to a situation different from that in which, in 
the real world, the residuals will actually correspond. 

The need for great flexibility in generating the IM becomes apparent: to account 
for the entire spectrum of possibilities occurring when the real analysis will be carried 
out. We therefore decided to build into the MATLAB routine the capability to vary 
randomly the initial phases (pk, the noise and the amplitude errors SAk- Concerning 
the error amplitudes of the harmonics, they can be randomly varied so that SAk spans 
[0, where (5A^°™ is the nominal mismodelled amplitude calculated taking into 

account Table V and Table IX; it means that we assume they are reliable estimates of the 
differences between the real data and the dynamical models of the orbital determination 
softwares, i.e. of the residuals. The MATLAB routine allows also the user to decide 
which harmonic is to be included in the IM; it is also possible to choose the length of the 
time series Tobs, the sampling step At and the amplitude of the noise. 

In order to assess quantitatively this feature we proceeded as follows. First, over a 
time span of 3.1 years, we fitted the IM with a straight line only, finding for a choice 
of random phases and noise which qualitatively reproduces the curve shown in ref. [66], 
the value of 38.25 mas for the root mean square of the post fit IM. The value quoted 
in ref.[66]is 43 mas. Second, as done in the cited paper, we fitted the complete IM with 
the LT plus a set of long-period signals (see subsect. 4) and then we subtracted the so 
adjusted harmonics from the original IM obtaining a "residual" simulated signal curve. 
The latter was subsequently fitted with a straight line only, finding a rms post fit of 12.3 
mas (it is nearly independent of the random phases and the noise) versus 13 mas quoted 
in [66]. A uniformly distributed noise with nonzero average and amplitude of 50 mas 
was used (see subsect. 4). These considerations suggest that the simulation procedure 
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adopted is reliable, replicates the real world satisfactorily and yields a good starting point 
for conducting our sensitivity analysis. 

4'4. Sensitivity analysis. - A preliminary analysis was carried out in order to evaluate 
the importance of the mismodelled solid tides on the recovery of /Ult- We calculated the 
average 



where [solid tides] denotes the analytical expressions of the mismodelled solid tidal per- 
turbations as given by eqs. (85)-(86). Subsequently, we compared it to the value of the 
gravitomagnetic trend for the same Tobs- The tests were repeated by varying At, Tobs 
and the initial phases. They have shown that the mismodelled part of the solid Earth 
tidal spectrum is entirely negligible with respect to the LT signal, falling always below 
1% of the gravitomagnetic shift of the combined residuals accumulated over the examined 



For the ocean tides and the other long-period signals the problem was addressed in a 
different way. First, for a given At and different time series lengths, we included in the 
IM the effects of LT, the noise and the solid tides only: subsequently we fitted it simply 
by means of a straight line. In a second step we have simultaneously added, both to the 
IM and the fitting model (FM hereafter), all the ocean tides, the solar radiation pres- 
sure SRP(657) signal and the odd zonal geopotential harmonic. We then compared the 
fitted values of /xlt and (5/xlt recovered from both cases in order to evaluate A/xlt and 
A^/xlt- The least squares fits [85, 86] were performed by means of the MATLAB rou- 
tine "nlleasqr" (sec, e.g., http://www.ill.fr/tas/matlab/doc/mfit.html); as (5/xlt 
we have assumed the square root of the diagonal covariance matrix element relative to 
the slope parameter. It simulates the formal experimental error which, in the context 
of the gravitomagnetic LAGEOS experiment, is far smaller than the systematic errors 
which, consequently, must be carefully assessed. Notice that the SRP(4,241) has never 
been included in the FM (See section 5), while Ki I = 3 p = 1 has been included for 
Tobs > 5 years. For both of the described scenarios we have taken the average for /jlt and 
f^/ULT over 1,500 runs performed by varying randomly the initial phases, the noise and 
the amplitudes of the mismodelled signals in order to account for all possible situations 
occurring in the real world, as pointed out in the previous section. The large number of 
repetitions was chosen in order to avoid that statistical fluctuations in the results could 
"leak" into A^Ilt and AJ/Ult and corrupt them at the level of 1%. With 1,500 runs the 
standard deviations on /Tlt and 5/iLT are of the order of 10""^ or less, so that wc can 
reliably use the results of such averages for our comparisons of 77LT(no tides) vs 77LT(all 
tides). 

Before implementing such strategy for different Tobsi At and noise we carehiUy ana- 
lyzed the At = 15 days, Tobs=4 years experiment considered in [15] by trying to obtain 
the quantitative features outlined there, so that we start from a firm and reliable basis. 
This goal was achieved by proceeding as outlined in the previous section and adopting a 
uniformly distributed random noise with an amplitude of 35 mas. 

Fig. 2 shows the results for A/Zlt obtained with At = 15 days and a uniform random 
noise with amplitudes of 50 mas and 35 mas corresponding to the characteristics of the 
real curves in [66, 15]. Notice that our estimates for the case Tobs = 4 years almost 
coincide with those in [15] who claim A/^tides < 4%. Notice that for Tobs > 7 years the 
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Influence of Eartfi tides on ^ 



noise 50 mas 
noise 35 mas 



5.5 6 6.5 

Obsen/ational period (years) 



Fig. 2. - Effects of the long period signals on the recovery of /ilt for At— 15 days and different 
choice of uniform randoni noise. Each point in the curves represents an average over 1,500 
runs performed by varying randomly the initial phases and the noise' s pattern. A/^lt is the 
difference between the least squares fitted value of /xlt when both the IM and the FM includes 
the LT plus all the harmonics and that obtained without any harmonic both in the IM and in 
the FM. 



effect of tidal perturbations errors falls around 2%. By choosing different At does not 
introduce appreciable modifications to the results presented here. This fact was tested 
by repeating the set of runs with At = 7 and 22 days. 

Up to this point we dealt with the entire set of long-period signals affecting the 
combined residuals; now we ask if it is possible to assess individually the influence of 
each tide on the recovery of LT. We shall focus our attention on the case At = 15 days, 
Tobs=4: years. 

In order to evaluate the influence of each tidal constituent on the recovery of /zlt 
two complementary approaches could be followed in principle. The first one consists of 
starting without any long-period signal either in the IM or in the FM, and subsequently 
adding to them one harmonic at a time, while neglecting all the others. In doing so it 
is implicitly assumed that each constituent is mutually uncorrclatcd with any other one 
present in the signal. In fact, the matter is quite different since if for the complete model 
case we consider the covariance and correlation matrices of the FM adjusted parameters 
it turns out that the LT is strongly correlated, at a level of | Corr{i, j)\ > 0.9, with 
certain harmonics, which happen to be mutually correlated too. These are 

• Ki, 1 = 2 (P=l,043.67 days; 1.39 cycles completed) 

• Ki, I = 2 (P=-569.21 days; 2.56 cycles completed) 

• Ki, 1 = 3, p = 2 (P=-336.28 days; 4.34 cycles completed) 

• K2, 1 = 3, p = 1 (P=-435.3 days; 3.35 cycles completed) 

• Solar Radiation Pressure, (P=657 days; 2.22 cycles completed) 

• perigee odd zonal C30, (P=821.79; 1.77 cycles completed) 

Their FM parameters are indeterminate in the sense that the values estimated are smaller 
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than the relative uncertainties assumed to be \JCov{i, i). On the contrary, there are other 
signals which are poorly correlated to the LT and whose reciprocal correlation too is very 
low and that arc well determined. These are 

• K2, 1 = 3, p = 2 (P=-211.4 days; 6.91 cycles completed) 

• 52, 1 = 3, p = 1 (P=-128.6 days; 11.3 cycles completed) 

• S2, I = 3, p = 2 (P=-97.9 days; 14.9 cycles completed) 

• Pi, I = 3, p = 1 (P=-166.2 days; 8.78 cycles completed) 

• Pi, 1 = 3, p = 2 (P=-118.35 days; 12.3 cycles completed) 

It is interesting to notice that over Tobs = 4 years the signals uncorrclatcd with the LT 
have in general described many complete cycles, contrary to these correlated with LT. 
This means that the signals that average out over Tobs decorrelate with LT, allowing the 
gravitomagnetic trend to emerge clearly against the background, almost not affecting the 
LT recovery. This feature has been tested as follows. In a first step, all the uncorrelated 
tides have been removed from both the IM (50 mas uniform random noise) and the 
FM, leaving only the correlated tides in. The runs were then repeated, all other things 
being equal, the following values were recorded: /Ilt = 1.2104, 5/iLT = 0.1006 with a 
variation with respect to the complete model case (/Ilt = 1-2073, (5/xlt = 0.1295 ) of: 
A/Tlt — 0.3%, AS^uT — 2.9%. Conversely, if all the signals with strong correlation are 
removed from the simulated data and from the FM, leaving only the uncorrelated ones in, 
we obtain: Jifjf = 1.1587, IJmjt = 0.0186 with Ajlfjf ~ 4.8%, A(5/xlt ^ 11%. Notice that 
the sum of both contributions for AjTuf yields exactly the overall value A/IiTf = 5.2% 
as obtained in the previous analysis (cf. fig. 2 for Tobs = 4 years). This highlights the 
importance of certain long period signals in affecting the recovery of LT and justify the 
choice of treating them simultaneously as it was done in obtaining fig. 2. Moreover, it 
clearly indicates that the efforts of the scientific community should be focused on the 
improvement of the knowledge of these tidal constituents. It has been tested that for 
^obs = 8 years all such "geometrical" correlations among the LT and the harmonics 
nearly disappear, as it would be intuitively expected. 



4'5. The effect of the very long periodic harmonics. - In this section we shall deal with 
those signals whose periods are longer than 4 years which could corrupt the recovery of 
LT resembling superimposed trends if their periods are considerably longer than the 
adopted time series length. 

In section 3 we have shown the existence of a very long periodic ocean tidal perturba- 
tion acting upon the perigee of LAGEOS II. It is the Ki I = 3 p = \ q = —1 constituent 
with period P ~ 1,851.9 days (5.07 years) and nominal amplitudes of —1, 136 mas. In 
ref. [87], for the effect of the direct solar radiation pressure on the perigee of LAGEOS 
II, it has been explicitly calculated, by neglecting the effects of the eclipses, a signal 
SRP(4,241) with P = 4, 241 days (11.6 years) and nominal amplitude of 6400 mas. The 
mismodelling on these harmonics, both of the form 5Asm(^t + (j)), amounts to 64.5 
mas for the tidal constituent, as shown in section 3.6, and to 32 mas for SRP(4,241), 
according to ref. [87]. 

About the actual presence of such semisecular harmonics in the spectrum of the real 
combined residuals, it must be pointed out that, over Tobs = 3.1 years [66], it is not 
possible that so low frequencies could be resolved by Fourier analysis. Indeed, according 
to refs.[88, 89], when a signal is sampled over a finite time interval Tobs it induces a 
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sampling also in the spectrum. The lowest frequency that can be resolved is 



(148) /„.i„ = 

i.e. a harmonic must describe, at least, half a cycle over Tobs m order to be detected in the 
spectrum, /min is called elementary frequency band and it also represents the minimum 
separation between two frequencies in order to be resolved. For our two signals we have 
f{Ki) = 5.39 X 10"^ cycles per day (cpd) and /(SRP) = 2.35 x 10"^ cpd; over 3.1 years 
/min = 4.41 X 10^'' cpd. This means that the two harmonics neither can be resolved 
as distinct nor they can be found in the spectrum at all. In order to resolve them we 
should wait for Tobs = 4.5 years which corresponds to their separation A/ = 3.04 x 10~^ 
cpd, according to eq. (148). In view of the potentially large aliasing effect of these two 
harmonics on the LT it was decided to include both Ki I = 3 p = 1 and SRP (4,241) in 
the simulated residual curve at the level of mismodelling claimed before. 

In order to obtain an upper bomid of their contribution to the systematic uncertainty 
in /iLT we proceeded as follows. First, we calculated the temporal average of the per- 
turbations induced on the combined residuals by the two harmonics over different Tobs 
according to 



(149) J = — / C26Asm{—t + (p)dt = C2 — (cos^ + sin </> sin r — cos cost), 

Jobs Jo ^ T 

with T = 27r2^. It turns out that, for given Tobs, the averages are periodic functions 
of the initial phases ^ with period 27r. Since the gravitomagnetic trend is positive and 
growing in time, we shall consider only the positive values of the temporal averages 
corresponding to those portions of sinusoid which are themselves positive. Moreover, 
notice that, in this case, one should consider if the perturbation' s arc is rising or falling 
over the considered Tobs: indeed, even though the corresponding averages could be equal, 
it is only in the first case that the sinusoid has an aliasing effect on the LT trend. The 
values of which maximize the averages were found numerically and for such values 
the maxima of the averages were calculated. Subsequently, these results in mas were 
compared to the amount of the predicted gravitomagnetic shift accumulated over the 
chosen Tobs by the combined residuals j/lt = 60.2 (mas/y)xTobs (y)- The results are 
shown in Table XIV and Table XV and, as previously outlined, represent a pessimistic 
estimate. 



Table XIV. - Effect of the averaged mismodelled harmonic K\ I = Z p = 1 on the Lense- 

Thirring trend for different Tobs- In order to obtain upper limits the maximum value for the 
average of the tidal constituent has been taken, while for the gravitomagnetic effect it has been 
simply taken the value yur x Tobs. 



Tobs (years) 


max C2 ('5'^^^) (mas) 


j/LT (mas) 


max C2 /2/LT (%) 


4 


5.6 


240.8 


2.3 


5 


0.3 


301 


0.09 


6 


3.3 


361.2 


0.9 


7 


4.8 


421.4 


1.1 
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Table XV. - Effect of the averaged mismodelled harmonic SRP(4, 241) on the Lense-Thirring 
trend for different Tobs • In order to obtain upper limits the maximum value for the average of 
the radiative harmonic has been taken, while for the gravitomagnetic effect it has been simply 
taken the value yur x Tpbs- 



Tobs (years) 


max C2 ("^t^") (mas) 


2/LT (mas) 


max C2 ("^t^") /2/LT (%) 


4 


9.1 


240.8 


3.7 


5 


8 


301 


2.6 


6 


6.8 


361.2 


1.8 


7 


5.6 


421.4 


1.3 



It is interesting to notice that an analysis over Tobs = 5 years, that should not be 
much more demanding than the already performed works, could cancel out the effect of 
the Ki I = 3 p = 1 tide; in this scenario the effect of SRP(4,241) should amount, at 
most to 2.6% of the LT effect. Notice also that for Tobs = 7 years, a time span sufficient 
for the LT to emerge on the background of most of the other tidal perturbations, the 
estimates of fig. 2 arc; compatible with those of Table XIV and Table XV which predict 
an upper contribution of 1.1% from the Ki I = 3 p = 1 and 1.3% from SRP(4,241) on 
the LT parameter /xlt- 

An approach, similar to that used in [90] in order to assess the influence of the 
eclipses and Earth penumbra effects on the perigee of LAGEOS II has been applied also 
to our case. It consists of fitting with a straight line only the mismodelled perturbation 
to be considered and, subsequently, comparing the slope of such fits to that due to 
gravitomagnetism which is equal to 1 in units of 60.2 mas/y. We have applied this method 
to Ki I = 3 p = 1 and SRP(4,241) with the already cited mismodelled amplitudes and 
by varying randomly the initial phases 4> within [— 27r, 2it]. The mean value of the fit' s 
slope, in units of 60.2 mas/y, over 1,500 runs is very close to zero. Concerning the upper 
limits of A/zlt derived from these runs, they agree with those released in Table XIV and 
Table XV up to 1-2 %. 

The method of temporal averages can be successfully applied also for the I = 2 m = 
18.6-year tide. It is a very long period zonal tidal perturbation which could potentially 
reveal itself as the most dangerous in aliasing the results for /xlt since its nominal am- 
plitude is very large and its period is much longer than the Tobs which could be adopted 
for real analysis. Ciufolini in ref.[14] claims that the combined residuals have the merit 
of cancelling out all the static and dynamical geopotential' s contributions of degree 
Z = 2,4 and order m = 0, so that the 18.6-year tide would not create problems. This 
topic was quantitatively addressed in a preliminary way. Indeed, in order to make com- 
parisons with other works, we have simply calculated for Tobs = 1 year the combined 
residuals with only the mismodelled amplitudes of the 18.6-year tidal perturbations on 
the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. This means that 
the dynamical pattern over the time span of such important tidal perturbation was not 
investigated. We did this by calculating the average over different Tobs- The results show 
clearly that the 18.6-year tide does not affect at all the estimation of /zlt if we adopt as 
observable the combined residuals proposed by Ciufolini. Indeed, for Tobs = 4 years, it 
turns out that the average effect will reach, at most, 0.08% of the gravitomagnetic shift 
over the same time span. 
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This feature of the 18.6-year tide is confirmed also by fitting with a straight hne only 
the sinusoid of this perturbation on the combined residual: the adjusted slope amounts, 
at most, to less than 1% of the gravitomagnetic effect for different Tobs- 

4'6. Conclusions. - The detection of the Lense-Thirring effect can be achieved by 
means of the cq. (5). Ciufolini in rcf. [14] claims that the static and dynamical terms of 
degree ? = 2, 4 and order m = of the geopotential do not affect his proposed combined 
residuals. 

In this section wc have tested this statement in regard not only to the Z = 2, 4 m = 
solid and ocean Earth tides, but also to the m = 1, 2 ones and to the other long-period 
perturbations which affect the combined residuals. 

Concerning the zonal tides, in order to compare our results to those released in [66], 
we have preliminarily used the nominal perturbative amplitudes by calculating eq. (5) 
on a time span of 1 year. As far as the I = 2 m = tides are concerned, the 18.6- 
year tide cancels out at a level of 10~^ only. But if we repeat the calculations with 
the mismodelled amplitudes over 1 year the accuracy of the cancellation grows to 10~^. 
This shows that the more accurate the dynamical models employed in building up the 
orbital residuals are, the more accurate the recovery of the Lense-Thirring becomes. The 
I = 4 zonal tides do not create problems. Also the Z = 3 m = tides, and this is an 
unpredicted feature, although at a lesser extent, cancel out at a level ranging from 10~^- 
10~^ if they are properly accounted in the residuals' construction. This points towards 
a better understanding, from both a theoretical and experimental point of view, of the 

I — 3, m = part of the tidal spectrum. The results presented in this section not only 
confirm the usefulness of the formula by Ciufolini for the / = 2, 4 m = tides, but also 
extend its validity to the Z = 3 m = part of the tidal response spectrum. 

In regard to the tesscral and sectorial tides, on one hand, by simulating the real 
residual curve in order to reproduce as closely as possible the results obtained for the 
At = 15 days Tobs = 4 years scenario published in [15] it has been possible to refine and 
detail the estimates for it. On the other hand, this procedure also extends them to longer 
observational periods in view of new, more sophisticated analysis to be completed in the 
near future based on real data analyzed with the orbit determination software GEODYN 

II in collaboration with the teams from the Joint Center for Earth Systems Technology 
at NASA Goddard Space Flight Center and at the University of Rome La Sapienza. 
Since such numerical analysis is very demanding in terms of both time employed and 
results analysis burden, it should be very useful to have a priori estimates which could 
better direct the work. This could be done, e.g., by identifying which tidal constituents 
/Ult is more or less sensitive to in order to seek improved dynamical models for use in 
GEODYN. 

As far as the perturbations generated by the solid Earth tides, the high level of 
accuracy with which they are known has yielded a contribution to the systematic errors 
in /iLT which falls well below 1%, so that they are of no concern at present. 

Concerning the ocean tidal perturbations and the other long-period harmonics, for 
those whose periods are shorter than 4 years, the role played by Tobs, At and the noise has 
been investigated. It turned out that At has no discernible effect on the adjusted value 
of HhT, while Tobs is very important and so is the noise. The main results for these are 
summarized in fig. 2, which tells us that the entire set of long-period signals, if properly 
accounted for in building up the residuals, affect the recovery of the Lense-Thirring effect 
at a level not worse than 4% — 5% for Tobs = 4 years; the error contribution diminishes 
to about 2% after 7 years of observations. 
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We have also shown which tides are strongly anticorrelated and correlated with the 
gravitomagnetic trend over 4 years of observations. The experimental and theoretical 
efforts should concentrate on improving these constituents in particular. This geometrical 
correlation tends to diminish as Tobs grows. This can be intuitively recognized by noting 
that the longer Tobs is, the larger the number of cycles these periodic signals are sampled 
over and cleaner the way in which the secular Lense-Thirring trend emerges against the 
background "noise". 

The ocean tide constituents Ki I = 3 p = 1 and the solar radiation pressure harmonic 

SRP(4,241) generate perturbations on the perigee of LAGEOS II with periods of 5.07 
years and 11.6 years respectively, so that they act on the Lense-Thirring effect as biases 
and corrupt its determination with the related mismodelling: indeed they may resemble 
trends if Tobs is shorter than their periods. They were included in the simulated residual 
curve and their effect was evaluated in different ways with respect to the other tides. 
An upper bound was calculated for their action and it turns out that they contribute to 
the systematic uncertainty in the recovery of /ilt at a level of less than 4% depending 
on Tobs and the initial phases i/i. The results are summarized in Table XIV and Table 
XV. An observational period of 5 years, which seems to be a reasonable choice in terms 
of time scale and computational burden, allows one to average out the effect of the 
Kil = 3p=l. 

The strategy followed for the latter harmonics has been extended also to the I = 
2 m = zonal 18.6-year tide. The obtained results confirm the claim in [14] that it does 
not affect the combined residuals. 

In conclusion, the strategy presented here could be used as follows. Starting from a 
simulated residual curve based on the state of art of the real analysis performed until 
now, it provides helpful indications in order to improve the force models of the orbit 
determination software as far as tidal perturbations are concerned and to perform new 
analysis with real residuals. Moreover, when real data will be collected for a given 
scenario it will be possible to use them in our software in order to adapt the simulation 
procedure to the new situation; e.g. it is expected that the noise level in the near future 
will diminish in view of improvements in laser ranging technology and modelling. Thus, 
we shall repeat our analysis for A/ittides when these new results become available. 

5. Measuring the gravitoelectric perigee advance with LAGEOS cind LA- 
GEOS II 

In order to test general relativity versus alternative metric theories of gravitation, fol- 
lowing a method introduced by Eddington (1922), Robertson (1962) and Schiff (1967), 
one can expand the spacctime metric (which describes the gravitational interaction) 
at the order beyond Newtonian theory (post-Newtonian) and then multiply each post- 
Newtonian term by a dimensionless parameter, the PPN parameters, to be experimen- 
tally determined. Using this method, Nordtvedt and Will have developed the PPN 
(Parameterized-Post-Newtonian) formalism [91], a powerful and useful tool for test- 
ing general relativity and alternative metric theories. Therefore, the PPN formalism 
is a post-Newtonian parameterized expansion of the metric tensor g and of the energy- 
momentum tensor T in terms of small known classical potentials. 

In the slow-motion and weak-field approximation, at the so-called post-Newtonian 
limit, the structure of known metric theories of gravity is identical apart from the nu- 
merical values of the PPN parameters which appear in the expansion of the metric 
coefficients [91, 8, 92]. Especially meaningful are the parameters (3 and 7, i.e the usual 
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Eddington-Robertson-Schiff parameters used to describe the "classical" tests of General 
Relativity. 7 measures how much space curvature is produced by unit rest mass and in 
the standard parameterized post-Newtonian gauge [3 accounts for the level of nonlinearity 
in the superposition law for gravity [32] . In general relativity (3 ~ 1 and 7=1. 

To-day a very accurate measurement of 7 was performed by using the Viking time 
delay [93]: the results was 7 = 1.000 ± 2 x lO^''. The quoted uncertainty allows for 
possible systematic errors. Other recent measurements exploit the gravitational bending 
of the electromagnetic waves at various wavelength [94] obtaining an uncertainty of the 
order of 10""^; e.g., according to the VLBI radio analysis by [95] 7 = 0.9996 ± 1.7 x 10~^, 
while in ref.[94] 7 = 0.997 ± 3 x 10~^ is quoted based on the astrometric observations 
of the electromagnetic waves deflection in the visible. An improvement in the accuracy 
of two orders of magnitude [94] is expected from the future GAIA astrometric mission 
[96] . Lunar laser ranging (LLR) measurements of the geodetic precession yields for 7 an 
accuracy of 10~^ [97]. However a careful discussion of the difii'erent terms in the lunar 
motion yields eventually 7 = 1.000 dz 5 x 10^"^ as the present LLR result [94]. 

The value of the parameter /3, which cannot be measured independently, can be 
obtained from the measured values of 7 and of some combinations of 7 and (3: r] = 
4/3 — 7 — 3 and are the most accurately determined. General relativity predicts 

?7 = and ^+^3^"^ = 1. For [3 in ref.[98] by using LLR data (5/3 = 4.7 x 10"'' is reported. 

The combination 77 measures the possible violation of the strong equivalence principle, 
i.e. the so-called Nordvedt effect [99]. 77 has been measured by analyzing the motion of 
the Moon via LLR data analyses [100, 97, 98]. The most recent determination of 77 by 
LLR is = 0.0002 ± 8 x 10"^ [98]. 

The combination ^"'"^3 ~^ is related to the well known pericenter shift of a test body 
induced by the Schwarzschild' s gravitoelectric part of the metric for a static, spherically 
symmetric distribution of mass-energy [32, 8]. This effect has been accurately measured 
for the Mercury perihelion shift in the field of the Sun by means of the echo delays of radar 
signals transmitted from Earth to Mercury [101] yielding ^^^3 ~^ = 1.005 ± 7 x 10~^. 
Inclusion of the probable contribution of systematic errors raises the uncertainty to about 
2x 10~^. In this case the major sources of systematic error are the poorly-known variation 
in topography on the planet's surface and the uncertainties in the radar scattering law 
[102, 103]. 

It should be considered that the interpretation of the perihelion advance of Mercury 

as a test of general relativity must cope with the uncertainty in the mass quadrupole 
moment of the Sun which also contributes to the perihelion advance. From the expres- 
sion of the perihelion precession and from the 1976 experimental uncertainties in w, one 
can easily calculate that for Mercury orbiting the Sun, any value of J20 larger than 
about 3x10"^' would disagree with the general relativistic prediction of 42.98 arcsec- 
onds/century; indeed, according to some authors measured values of J20 may be as large 
as ~ 5.5xl0~^ (for a comprehensive discussion see refs. [8, 104]). Laboratory experiments 
have been also proposed [105] to measure the solar mode / = 1, to got information on the 
rotation of the core of the Sun from the I = 1 rotational frequency splitting. Nevertheless, 
according to more recent determinations of J20 using the analysis of the Sun's pressure 
modes, both from the ground-based network of observatories and the space based SOHO, 
J20 = (2.3±0.1)xl0-'^ [106]; in ref. [107] a theoretical estimate of J20 = (2.0±0.4) x lO"'^ 
is reported. Thus, according to these values of J20 the observed perihelion advance of 
Mercury is well in agreement with the general relativistic predictions. 

In order to measure J20, among other astrodynamical and relativistic parameters. 
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several space missions have been proposed using a variety of techniques: the most recent 
and promising are SORT, IPLR and ASTROD [108]. In [109] a measurement of the 
LAGEOS laser ranged satellite perigee shift in the field of the Earth is quoted, but the 
accuracy amounts only to 2 x 10~^. More accurate measurements of ^"'"^3 ~^ might be 
performed in the future by means of an ESA Mercury orbiter [110, 111]. 

In this section wc explore the possibility of performing a complementary measurement 
of ^"""^g ~^ in the gravitational field of the Earth by using some suitable combinations of 
the orbital residuals of the presently (or proposed) existing spherical geodetic laser ranged 
satellites with particular emphasis to LAGEOS and LAGEOS II in order to exploit the 
relevant experience obtained with the gravitomagnetic measurements in [112, 66, 15]. 
Moreover, as a secondary outcome of the proposed experiment, by combining the accurate 
determinations of the Nordtvedt parameter rj and of ^"'"^3 ~^ from SLR it would be 
possible to obtain independent values for 7 and /3. 

This section is organized as follows. In subsect. 1 we compare the present experi- 
mental accuracy in satellite laser-ranging measurements with the rclativistic expressions 
of the perigee shift of the two LAGEOS satellites. In subsect. 2 we analyze some of 
the most important sources of systematic errors: the even zonal harmonics of the static 
part of the geopotential, the Earth solid and ocean tidal perturbations, the direct solar 
radiation pressure and the mismodelling of the inclination. Subsect. 3 is devoted to the 
discussion of the obtained results and to the conclusions. 

5'1. The relativistic perigee precession of LAGEOS type satellites. - As known, in 
the slow-motion and weak-field approximation, the Schwarzschild metric generated by 
a static, spherically symmetric distribution of mass-energy induces an additional post- 
Newtonian "gravitoelectric" force which acts on the orbit of a test body by shifting its 
pericenter; in the PPN formalism the pericenter rate can be written as 

• 3nGM 2 + 27-/3 
= cMl - e^) ^ • 

In the following we define v = ^"'"^3 ~^ . General Relativity predicts that the perigee 
shifts for LAGEOS and LAGEOS II amount to 3,312.35 mas/y and 3,387.46 mas/y, 
respectively. 

Following ref.[14] the actual experimental precision allows for detecting such rates 
for both LAGEOS. Indeed, for the perigee the observable quantity is r = eauj and at 
present its measurement error amounts to about Sr^xp < 1 cm for the two LAGEOS, 
over several orbits and for a given set of force models (i.e. not including modelling 
errors). Since the LAGEOS eccentricity is ej = 4.5 x lO^'' the accuracy in detecting 
the perigee is Suil^p — Sve^p/ (eiai) ~ 37 mas. So, over 1 year the relative accuracy of 
the measurement of the relativistic perigee shift would be ~ 1 x 10~^. For LAGEOS 
II this measurement accuracy is better than that of LAGEOS indeed the LAGEOS II 
eccentricity is en = 1.4 x 10^^, so that ^w^xp — 12 mas; this may yield an accuracy 
of about 3 X 10~^ over 1 year. These considerations rule out the possibility of directly 
using the perigee of the other existing spherical geodetic laser-ranged satellites Etalon- 
1, Etalon-2, Ajisai, Stella, Westpac-1 because their eccentricities are even smaller than 
that of LAGEOS. On the contrary, Starlette has an eccentricity of the order of 2 x 10~^; 
however, since it orbits at a lower altitude it is more sensitive than the LAGEOS satellites 
to atmospheric drag and to Earth's zonal harmonics, so that it would be difficult to 
process its data at an acceptable level of accuracy. Accordingly, we will focus on the 
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perigee of LAGEOS and especially of LAGEOS II in order to accurately detect the 
gravitoelectric relativistic shift in the gravitational field of Earth. 

5'2. The systematic errors. - The most important source of systematic error in such 
measurements is represented by the mismodelling induced by the even zonal harmonics 
of the Earth gravitational field [1] on the classical perigee precession. By using the 
covariance matrix of the EGM96 Earth gravity model [19] and adding the correlated terms 
in a root-sum-square fashion up to degree Z = 20 we obtain for LAGEOS a systematic 
error Si'/i'zomds = 8.1 x 10~^, whereas for LAGEOS II we have (5i^/i^zonais = 1-5 x 10~^. 

Since the major source of uncertainty lies in the first two even zonal harmonics 6J2 and 
SJ4, following [14] for the Lense-Thirring LAGEOS experiment, we search for suitable 
combinations of orbital residuals of the existing SLR satellites in order to eliminate most 
static and dynamical even zonal terms of the geopotential. In Table XVI we report the 
most promising combinations: their general form is 

N 

(151) 6d}^^ + ^ CjdCl^ + cn+i^w' = xgr V, 

J=l 

in which N is the number of the nodes of different SLR satellites employed, xqr is the 
slope, in mas/y, of the relativistic trend to be measured and 5v / v.^on&Xs is the systematic 
error induced by the even zonal harmonics up to degree / = 20 calculated with EGM96 
covariance matrix. It is important to stress that the use of LAGEOS, due to their altitude, 
makes our measurement substantially insensitive to the errors in the zonal harmonics of 
degree I > 20, so that our estimates of Sv/vzonais presented here are valid even in the 
case that the EGM96 covariance matrix for higher degrees I > 20 would not be accurate 
enough. In addition to LAGEOS and LAGEOS II, we have only considered Ajisai since 



Table XVI. PPN combined residuals. 











J 






K 


Cl 


C2 


C3 


C4 


XGR (mas/y) 


^t'/l'zonalB 


1 


-0.868 


-2.855 








3,387.46 


6.59 X 10"^ 


2 


-2.514 


-4.372 





2.511 


11,704.92 


1.1 X 10"^ 


3 


-1.962 


-3.693 


0.0366 


1.370 


7,928.51 


8.1 X 10"* 



it is well tracked, contrary to, e.g., the Etalons, and it would be less demanding than 
the other satellites to reduce its laser ranged data to a level of accuracy comparable 
to that of LAGEOS and LAGEOS II. Moreover, the other SLR satellites orbit at lower 
altitudes, therefore they are more sensitive to the higher degree terms of the geopotential. 
Consequently, as confirmed by numerical calculations, the inclusion in their data in the 
combined residuals would increase Sv/vzonais,- 

Note that the systematic error induced by the mismodelled secular rates of the even 
zonal harmonics of the geopotential is really critical because it can be considered as an 
unavoidable part of the total error in the experiment. Indeed, the resulting aliasing trend 
cannot be removed from the data and nothing can be done about it apart from assessing 
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as more reliably as possible the related error. (^) 

It is important to point out that the values quoted in Table XVI for Sv / VzonaU will 
be reduced in the near future when the data from the CHAMP mission will be released. 

In regard to the evaluation of the impact of other sources of systematic errors in 
this section we consider in detail only the combination 1 of Table XVI so to exploit 
the background acquired with the Lense-Thirring LAGEOS experiment [112, 66, 15]. 
Moreover, reducing the Ajisai' s data to an acceptable level of accuracy for our measure- 
ment would neither be a trivial nor an immediate task to be performed and the inclusion 
of the perigee of LAGEOS would raise the experimental error. E.g., for combination 
2 of Table XVI, over 1 year, the impact of the error in measuring the perigee rate of 
LAGEOS amounts to 2.5 x (3.1 x 10"^) = 7.7 x 10"^ while for combination 3 it is 
1.370 X (4.7 X 10-3) = 6.4 x 10"^. 

5'3. Discussion and conclusions. - In Table XVII we summarize the results obtained 
for a 8 years long time span with 7 days arc lengths. The results and the methods of 
sections 3 and 4 and of ref.[87, 113] have been used. 

Table XVII. - Preliminary error budget: Tabs = 8 years, At = 7 days. 



Even zonal harmonics 
Js geopotential 
Tides 

Non-gravitational effects 
Measurement error in LAGEOS II perigee 

Total systematic error 7.3 x 10 



6.59 X 10" 
3.2 X 10"* 
4.4 X lO""* 
1 X 10"* 
4 X 10"* 



In regard to the non-gravitational errors, according to the results of ref. [87] for the 
direct solar radiation pressure and the Earth's albedo and to the results of ref. [113] 
for the thermal thrust perturbations and the asymmetric reflectivity, the corresponding 
uncertainty would amount to almost 1 x 10"^ over 7 years. In obtaining this result a very 
pessimistic approach has been adopted by assuming a mismodeling of 20% for all the 
perturbing effects except for the direct solar radiation pressure which has been assumed 
to have an uncertainty at the 0.5% level. However, we stress that only the Earth's thermal 
thrust, or Yarkovski-Rubincam effect, induces a mismodeled linear trend whose impact 
would amount to about 1 x 10~^: the other non-gravitational forces are time-varying 
with known periodicities and can be fitted and removed from the signal as done for the 
tidal perturbations. 

In assessing the total systematic error we have accounted for the fact that the gravi- 
tational errors are not independent simply by summing them up. Then, we have added 
them and the other independent sources of errors in a root-sum-square fashion. Note 
that the error in the LAGEOS II perigee refers to 1 year only by assuming an error of 1 
cm in its radial position. Moreover, the estimate of the statistical formal error has been 



(^) About combination 3, it should be noted that in order to obtain more rehable and accurate 

estimates of the systematic error due to the even zonal harmonics of the geopotential it should 
be better to extend the calculations to higher degrees than I = 20 due to the sensitivity of Ajisai 
to such higher degree terms. 
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obtained by considering the noise level in the LAGEOS data of 1998 Lense-Thirring ex- 
periment; in the near future and for longer time spans these uncertainties will be reduced. 
Moreover, it is important to emphasize that the data from the CHAMP and GRACE 
missions will soon yield a notable reduction of the systematic error due to the higher 
degree static even zonal harmonics of the geopotential and of the solid and ocean tides 
systematic errors as well. In particular, a very significant reduction of the systematic 
errors will take place due to the static part of the geopotential, which has turned out to 
be the most important source of uncertainty. 

The results obtained for the combination 1 examined here together with those released 
in ref.[98] for the combination r] would allow to measure (3 = jT] + ji/ + j independently 
of other measurements of 7 with 3.3 x 10~^ accuracy over 8 years; this result, which is 
of the same order of magnitude of that obtained with the radar ranging technique [102], 
should be compared with the most recent dp = 4.7 x 10^^ obtained from the LLR data 
[98]. The parameter 7, given by 7 = jr]+ — |, could be measured less precisely over 
the same time span: ^7 1.2 x 10~^. 

6. — The LARES mission revisited: some alternative scenarios 

6T. Introduction. In order to measure the Lense-Thirring drag with an accuracy of 
few percent, in [12] it was proposed to launch a passive geodetic laser-ranged satellite- 
the former LAGEOS III - with the same orbital parameters of LAGEOS apart from its 
inclination which should be supplementary to that of LAGEOS. 

This orbital configuration would be able to cancel out exactly the classical nodal 
precessions provided that the observable to be adopted is the sum of the residuals of the 
nodal precessions of LAGEOS III and LAGEOS 

(152) 6tl™ + = 62/ZLT. 

Later on the concept of the mission slightly changed. The area-to-mass ratio of LA- 
GEOS III was reduced in order to make less relevant the impact of the non-gravitational 
perturbations and the eccentricity was enhanced in order to be able to perform other 
general relativistic tests: the LARES was born. The most recent error budget of the 
LARES Lense-Thirring experiment claim an accuracy of the order of 3%. Unfortunately, 
at present we do not know if the LARES mission will be approved by any space agency. 

In this section we investigate the possibility of further reducing the error in this impor- 
tant proposed mission. This section is organized as follows. In subscct. 2 we analyze in 
detail the impact of the unavoidable injection errors in the orbital parameters of LARES 
on the systematic error induced by the mismodelling in the even zonal harmonics of the 
geopotential according to the most recently released Earth gravity model. Moreover, 
in subsect. 3 we propose an alternative configuration which should be able to reduce 
this error by one order of magnitude. It adopts as observable a suitable combination 
of the orbital residuals of the nodes of LAGEOS, LAGEOS II and LARES, and the 
perigees of LAGEOS II and LARES. It presents also the important advantage that it 
is almost insensitive to the errors in the inclination of LARES, contrary to the original 
LAGEOS/LARES only configuration. The negative implications of placing the LARES 
in a low-altitude polar orbits are examined in subsect. 4. Subsect. 5 is devoted to the 
conclusions. 
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6'2. The impact of the even zonal harmonics of the geopotential on the original LARES 
mission. - Let us calculate the systematic error induced by the mismodelling in the 
even degree zonal coefficients J2, J4,... of the geopotential on the sum of the classical 
precessions of the nodes of LAGEOS and LARES. It is important to stress that it is the 
major source of systematic error and it cannot be eliminated in any way. We will use 
the covariance matrix of the Earth gravity field model EGM96 [19] by summing up in a 
root sum square fashion the correlated contributes up to degree I = 20(^°). The relative 
error obtained by using the nominal values amounts to 

(153) ^ = 3 X 10-3. 

MLT zonals 

It is not equal to zero because wc have assumed eLR = 0.04 while glageos = 0.0045. 
If it was eLR = clageoSj then the classical nodal precessions would be exactly equal 
in value and opposite in sign and would cancel out. Notice that the coefficients with 
which SCl™ and 60^ enter the combination of eq. (152) do not depend on any orbital 
parameters: they are constants equal to 1. Moreover, eq. (152) is affected by all the 
classical nodal precessions, including those induced by J2 and J4 which, instead, are 
cancelled out a priori in the combination used in the LAGEOS experiment [14]. They 
are the most effective in aliasing the Lense-Thirring precessional rates. 

Now we will focus on the sensitivity of to the unavoidable orbital injection 

A^LT zonals 

errors in the possible orbital parameters of LARES. For a former analysis see ref.[114]. It 
is particularly interesting to consider the impact of the errors in the inclination and the 
semimajor axis. The ranges of variation for them have been chosen in a very conservative 
way in order to take into account low-precision and low-costs injection scenarios. 

From fig. 3 it is interesting to notice that the minimum value of the systematic zonal 
error, which amounts to 2.1 x 10"^, does not correspond to ^lr = 70 deg but it is obtained 
for a slightly smaller value. It is possible to show that for rlr = eLACEOS the minimum 
is and that it is attained at ilr = 70 deg. The maximum error amounts to 1.6 x 10"^. 
This suggests that the original LARES project is rather sensitive to small departures 
of iLR from its nominal value. Fig. 4 shows that even more relevant is the sensitivity 
to the LARES semimajor axis. Also in this case the minimum is attained at a value 
of ttLR smaller than the nominal clr = 12, 270 km. For glr = eLAGEOS the minimum 
error amounts to and it is obtained for olr = 12, 270 km, as expected. In obtaining 
fig. 4 we have accounted for the dependence of the nodal Lense-Thirring precession on 
a by varying, accordingly, the slope of the general relativistic trend. The sensitivity to 
eccentricity variations is less relevant: indeed, by varying it from 0.03 to 0.05 the relative 
systematic zonal error passes from 1.6 x 10"^ to 4.6 x 10"^ 

6'3. An alternative LARES scenario. - Here we will look for an alternative observable 
involving the orbital elements of LARES which 

• yields a smaller value for the systematic error due to the mismodclled even zonal har- 
monics of the geopotential than that of the simple sum of the nodes of LAGEOS and 
LARES 



i}'^) It is important to notice that using satellites of the LAGEOS family allows one to obtain 
reliable estimates with EGM96. Indeed, the LAGEOS' orbits are not affected by the terms of 
degree I > 20, for which the EGM96 covariance matrix elements are determined less accurately. 
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SnU 5n'"'^=62ji ■ sensitivity to orbital injection errors in inclination 



constant combination coefficients = 1 ; a, „=12,270 km; e, „=0.04 




69 69.2 69.4 69.6 69.8 70 70.2 70.4 70.6 70. f 

LARES inclination (deg) 



Fig. 3. - Influence of the injection errors in the LARES inclination on the zonals' error. 



5n + S£i =X|i ■ sensitivity to orbital injection errors in LARES semimajor axis 
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Slope of the gravitomagnetic trend: sensitivity to orbital injection errors in LARES serr 
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Fig. 4. - Influence of the injection errors in the LARES semimajor axis on the zonals' error. 



• is less sensitive to the departures of the possible orbital elements of LARES from the 
nominal values. The first requirement could be implemented by setting up a suitable 
orbital combination which cancels out the contributions of as many mismodelled even 
zonal harmonics as possible, following the strategy of the LAGEOS experiment outlined 
in [14]. To this aim we will consider only the satellites of the LAGEOS family both 
because they are the best laser-ranged tracked targets and because the gravitational and 
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non-gravitational perturbations affecting their orbits have been extensively and thor- 
oughly analyzed. 
Our result is 

with 



(155) Cl = 6 X 10-^ 

(156) C2 = 9.83 X 10-\ 

(157) C3 = -1 X 10-^ 

(158) C4 = -2 X 10"^ 



It is important to notice that the coefficients given by cqs. (155)-(158) depend on the or- 
bital parameters of the satellites entering the combination and, among them, of LARES. 
The values released here are calculated for the nominal LARES parameters, as is the 
case for the slope in mas/y of the general rclativistic trend. The observable of eq. (154) 
allows one to cancel out the static and dynamical contributions of the first four even 
zonal harmonics. The relative systematic error due to the J2„, n > 5, according to 
EGM96 up to degree I = 20, amounts to: 

(159) ^ = 2 X 10-^ 

MLT zonals 

which is one order of magnitude better than the result of eq. (153). 

Fig. 5 and fig. 6 show the important achievements realized in reducing the sensitivity 
of the proposed combined residuals to the orbital injection errors in the LARES orbital 
elements. In obtaining fig. 5 and fig. 6 we have accounted for the dependence on olr 
and iLR of both the coefficients and the Lcnsc-Thirring precessions: it turns out that the 
variations in the slope of the general rclativistic trend arc very smooth with respect to 
the nominal value of 61.8 mas/y amounting to few mas/y. Now the values of the zonals' 
error are much more close to the nominal one given by eq. (159). Also in this case, 
the minima are attained at slightly different values of the LARES orbital elements with 
respect to the nominal ones. It is interesting to notice in fig. 5 that over a 3% variation of 
Glares the error due to the mismodeled zonal harmonics remain almost constant, while 
over a 5% variation of aLARES it changes of 1 order of magnitude, as it turns out from 
fig. 6. It is worthwhile noting that the time-varying gravitational and non-gravitational 
orbital perturbations which affect the proposed combined residuals are depressed by the 
small values of the coefficients with which some orbital elements enter the combination. 
Moreover, the observational error in the LAGEOS II and LARES perigees, which are 
difficult to measure for low eccentric satellites as LAGEOSdue to the small value of their 
eccentricity, would have an impact of the order of 1 x 10~* by assuming an uncertainty 
of 1 cm over 1 year in the satellite's position. 

6'4. Conclusions. - If analyzed from the point of view of the impact of the systematic 
error induced by the mismodelling in the even zonal harmonics of the geopotential, which 
is the most important source of systematic error, the original LARES mission seems to 
be affected by a certain sensitivity to the unavoidable departures of the LARES orbital 
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Fig. 5. - Alternative combined residuals: influence of the injection errors in the LARES incli- 
nation on the zonals' error. 
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Fig. 6. - Alternative combined residuals: influence of the injection errors in the LARES semi- 
major axis on the zonals' error. 



parameters from their nominal values due to orbital injection errors. In such a refined 
experiment, which would compete with the ambitious Stanford GP-B mission and its 
claimed global 1% accuracy level, it could be a serious drawback. It could be enhanced 
if LARES will be finally put in orbit with a low-cost launcher which, inevitably, would 
induce relatively large injection errors. 
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The adoption of the alternative combined residuals proposed here, including also 
the node of LAGEOS II and the perigees of LAGEOS II and LARES, would reduce 
by one order of magnitude the systematic error due to the even zonal harmonics of 
the geopotential passing from 0.3% to 0.02%, and would reduce greatly the sensitivity 
of such result to errors in the LARES orbital parameters. This would yield to less 
stringent requirements on the quality and the costs of the launcher to be adopted. It 
is very important to notice that the estimates presented here are based on the most 
recent Earth gravity model EGM96. When the new data on the terrestrial gravitational 
field from the CHAMP and GRACE missions will be available, the systematic zonal 
error will greatly reduce. The impact of the errors related to the quality of laser data 
will further reduce in the near future as well. Preliminary estimates of the standard 
statistical error in the solve- for least square parameter /ilt, based on the present models 
of the time-dependent LAGEOS perturbations worked out in section 3 and the noise level 
reported in the Lense-Thirring LAGEOS experiment, yield a value of the order of 10~^. 
If confirmed by further analysis of the impact of the time-dependent gravitational and 
non-gravitational perturbations, this result is very important because the total error of 
this improved version of the LARES mission would be < 1%. 

7. — Conclusions and recommendations for future work 

The general rclativistic features considered here are 

• The gravitomagnetic LAGEOS experiment devoted to the measurement of the 
Lense-Thirring drag of inertial frames, currently carried on by analyzing the combined 
residuals of the nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II. The 
assessment of the systematic error induced by the time-dependent part of the terrestrial 
gravitational field has been extensively discussed so to extend its evaluation also to future 
measurements over longer time spans. The role of the proposed LARES in a new context 
is also taken into account. The proposed new observable would allow to reduce by one 
order of magnitude the total systematic error 

• The pericenter shift of a test particle due to the Schwarzschild's gravitoelectric part 
of the gravitational field of a central, static body and the possibility of measuring it at 
a 10"'' level in the field of the Earth with LAGEOS, LAGEOS II, and, perhaps, other 
SLR satellites 

It is very important to stress that in the near future, when the new and more precise 
data relative to the static and dynamical parts of the Earth's gravitational field from the 
CHAMP and GRACE missions will be available, it will be possible to reduce sensibly 
the uncertainties of many source of systematic errors of gravitational origin. This is 
particularly true for the even zonal harmonics of the static part of the geopotential. 
In regard to the gravitomagnetic LAGEOS experiment, this would probably allow to 
reduce the total geopotential error from the present 12.9% to some percent. Moreover, 
the zonals error in the PPN measurement could reach the level of 10"''— 10~^ according to 
the observable adopted. These improvements would extend also to the LARES project. 
So, it could be possible to proceed in the future by updating the orbit determination 
softwares like GEODYN II with the new data of the geopotential and, subsequently, by 
reanalyzing the data of the gravitomagnetic LAGEOS experiment. In the framework 
of the gravitomagnetic LAGEOS experiment, it could be useful to look for alternative 
combinations of orbital residuals involving also other SLR satellites in order to reduce, 
on one hand, the systematic error due to the geopotential, on the other, the systematic 
error due to the non-gravitational perturbations. It would be the follow-on of the work 
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presented in ref.[114]. A first attempt can be found in ref.[115]. 
Appendix A. 

The classical orbital precessions 

Here we show the explicitly calculated coefficients 

^^class 



(A.l) O. 



2n = 



d{J2n) 

and 

(A.2) uj,2„ 



d{J2n) 

of the satellites' classical nodal and apsidal precessions due to the even {I = 2n, n = 
1,2,3...) zonal (m = 0) harmonics of the geopotential up to Z = 20. Recall that Ji = 

— C;o, I = "2,71, n = 1,2, 3... where the unnormalizcd adimcnsional Stokes coefficients 
of degree I and order m can be obtained from the normalized Cim with 

(A.3) Cim = Ni^Clm- 

In it 

' (21 + 1){2 - Som){l - m)\ 



(A.4) Ni^ 



(l + m)] 



For the general expression of the classical rates of the near Earth satellites' Keplerian 
orbital elements due to the geopotential ciciass, edass, «ciass, S^ciass, Wciass, class, see 
ref. [1]. The coefficients Q.2n and Lu.2n arc of crucial importance in the evaluation of 
the systematic error due to the mismodelled even zonal harmonics of the geopotential; 
moreover, they enter the combined residuals' coefficients Ci, i = 1,2...N of sect. 5 and 
sect. 6. See also ref. [14]. Since the general relativistic effects investigated are secular 
perturbations, we have considered only the perturbations averaged over one satellite' s 
orbital period. This has been accomplished with the condition I — 2p + q = 0. Since 
the eccentricity functions Gipq are proportional to e' for a given value of / we have 
considered only those values of p which fuffil the condition l — 2p + q = with q = 0, i.e. 
p = |. This implies that in the summations 

and 

; +OC , „ 

(A.6) ^F,, y: ^ 

p=0 g=— oo 
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involved in the expressions of the classical rates we have considered only -F;qi and Gjig. 
Moreover, in working out the GjiQ we have neglected the terms of order 0{e^) with 
k>2. 

The nodal coefficients, proportional to 
(A.7) ^ g G,„Y^'^ 

sin 5 '^^ 



sin i ^ di 
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The coefficients of the classical perigee precession are much more involved because 
they are proportional to 
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The second set is given by 
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The results previously obtained can be used in working out explicitly the contributions 
of the mismodelled classical nodal and apsidal precessions up to degree Z = 20 of the 
existing spherical passive laser-ranged geodetic satellites and of LARES. They are of the 
form (5f2(2n) = ^.2n X SJ2n, n = 1,2, ...10 and 5uj^2n) = ^'.2n X 5J2m n = 1,2, ...10. The 
values employed for SJ2n = V^rT+T x 5C2n 0, n = 1,2, ...10 would be those quoted in, 
e.g., EGM96 model. 



Appendix B. 

Some useful parameters used in the text 



The data employed for the terrestrial space environment and the LAGEOS satellites 
are in the following table. In it s is the angle between the ecliptic and the equatorial 
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plane, is the Earth mean equatorial radius, G is the Newtonian gravitational constant, 
GMq and 6{GM^) are the Newtonian gravitational constant times the Earth's mass and 
its error according to lERS standard [56], J2, J4, (5J2 and SJ^ are the first two even zonal 
geopotential's coefficients and their errors according to EGM96 [19], J© is the Earth's 
angular momentum, w© is the mean Earth angular velocity. 

The conversion factor from rad/s to mas/y is 1 rad/s=6.509 x 10^^ mas/y. 



Table XVIII. - Earth's and LAGEOS parameters used in the text. 



Parameter 


Numerical value 


Units 


e 


23.44 


deg 


Re 


6, 378 X 10® 


cm 


G 


6.67259 X 10"^ 


312 

cm' g s 


GM^ 


3.986 X 10^° 


cm"* s 


S{GM<s) 


8 X 10" 


cm"* s 


J2 


1.0826 X 10"^ 




J4 


-1.6194 X 10"® 




5J2 


7.9626 X 10"" 




SJ4 


3.126 X 10"^° 






5.9 X 10"" 


g cm^ s"^ 




3.3 X 10^ 


cm 


UJQ 


7.29 X 10"® 


rad s"^ 


G/c^ 


7.42 X 10"^^ 


cm g"'^ 


{GJe)/c' 


4.37 X 10^^ 


cm^ s"^ 


{GM^)/^ 


4.43 X 10"^ 


cm 


aLAGEOS 


1.2270 X 10^ 


cm 


dLAGEOS II 


1.2163 X 10^ 


cm 


OlARES 


1.2270 X 10^ 


cm 


eLAGEOS 


0.0045 




eLAGBOS II 


0.014 




eLARES 


0.04 




ILAGEGS 


110 


deg 


ILAGBOS II 


52.65 


deg 


jlares 


70 


deg 


nLAGEOS 


4.643 X 10"" 


s"i 


WLAGEOS II 


4.710 X 10"" 


s"i 


riLARES 


4.643 X 10"" 


s"i 



Appendix C. 
WEB resources 

• http: / /wugrav. wustl.edu/People/CLIFF/tegp. html, the experimental grav- 
ity web page. 

• http: / /www. livingreviews.org/Articles/Volume4/2001-4will/index. html, 

recent review of the current status of the empirical basis of General Relativity. 



68 L. lORIO 

• http://www.laeff.esa.es/eng/laefF/activity/lageos3.html, on LAGEOS III 
project. 

• http://w3.ing.uniromal.it/tildespacedpt/lares.html, on LARES project. 

• http: / /earth. agu.org/ revgeophys/mgirshaOl /nodel.html, on the force mod- 
els acting on LAGEOS satellites. 

• http://einstein.stanford.edu/index.html, on the GP-B mission. 

• http://www.nas.edu/ssb/gpbexe.html, on the GP-B mission. 

• http://op.gfz-potsdam.de/champ/, on the CHAMP mission. 

• http://op.gfz-potsdam.de/grace/, on the GRACE missions 

• http://ilrs.gsfc.nasa.gov/, the International Laser Ranging Service web site. 

• http://www.ee.surrey.ac.uk/SSC/SSHP/, small satellites web site. 

• http: / /cddisa.gsfc.nasa.gov/926/egm96/egm96.html, the EGM96 Earth grav- 
ity model. 

• http://co-ops.nos.noaa.gov/restlesl.html, on the phenomenon of Earth tides. 

• http:/ /www. astro. oma.be/Dl/EARTH_TIDES/wgtide. html, home page of 
the Working Group of Theoretical Tidal Model of the Earth Tide Commission. 

• http://www.ill.fr/tas/matlab/doc/mfit.html, useful collection of MATLAB 
least squares programs and other routines. 
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